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ABSTRACT 

Multispectral  Passive  Microwave  Correlations  with  an  Antecedent 
Precipitation  Index  using  the  Nimbus  7  SMMR.  (August  1984) 
Gregory  Delfin  Wilke,  B.G.S.,  University  of  Nebraska  Omaha 
M.S.,  University  of  Nebraska  Omaha 
Chairman  of  Advisory  Committee:  John  F.  Griffiths 

Analysis  of  the  passive  microwave  brightness  temperatures  from 
the  Scanning  Multichannel  Microwave  Radiometer  (SMMR)  aboard  Nimbus  7 
can  infer  soil  moisture  through  an  Antecedent  Precipitation  Index 
(API).  This  investigation  involves  correlation  analysis  between  five 
passive  microwave  brightness  temperature  wavelengths  (0.81,  1.36, 

1.66,  2.80,  and  4.54  cm)  in  two  polarizations  spanning  a  one  year 
period.  Data  analysis  was  done  by  a  temporal  and  areal  approach,  for 
both  the  normalized  brightness  temperatures  and  the  differences  in  the 
polarized  channels,  for  all  five  SMMR  wavelengths.  Three  transforma¬ 
tions  of  the  polarized  brightness  temperatures  were  used  to  evaluate 
the  ability  of  the  microwave  data  to  infer  soil  moisture  through 
correlations  with  the  API.  These  transformations  included  the 
horizontally  polarized  normalized  brightness  temperature,  the  dif¬ 
ference  between  the  polarized  channels,  and  the  normalized  difference 
between  the  polarized  channels.  It  was  evident  that  all  three  indices 
had  very  high  correlations  with  API,  however  the  correlation  coef¬ 
ficient  for  the  normalized  brightness  temperature  data  consistently 
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was  slightly  higher  than  the  values  for  the  other  two  indices.  The 
conclusion  is  that  polarized  data  can  infer  soil  moisture  through  the 
API  without  the  use  of  surface  data. 
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CHAPTER  I 

INTRODUCTION 

Purpose 

The  purpose  of  the  research  is  to  determine  the  relationship 
between  multispectral  dual  polarized  satellite  passive  microwave 
brightness  temperatures  and  the  Antecedent  Precipitation  Index  (API) 
over  uniform  agricultural  areas  in  the  Southern  Great  Plains  of  the 
United  States.  The  API  is  calculated  from  actual  precipitation  obser¬ 
vations  and  simulates  soil  moisture,  therefore,  soil  moisture  content 
can  be  inferred  through  knowledge  of  the  API  (Blanchard  et_^K,  1980). 
The  specific  instrument  used  to  obtain  the  data  was  the  Scanning  Multi¬ 
channel  Microwave  Radiometer  (SMMR)  aboard  the  Nimbus  7  satellite. 
Specific  objectives  included  the  investigation  of  three  transformations 
of  both  the  vertical  and  horizontal  brightness  temperatures  to  deter¬ 
mine  the  applicability  of  both  normalized  brightness  temperature  and 
polarized  channel  data  to  infer  soil  moisture. 

Need  for  Soil  Moisture  Information 

Knowledge  of  large-scale  patterns  of  soil  moisture  is  needed  to 
evaluate  the  condition  of  the  soil  or  crop  and  is  essential  to  the 
disciplines  of  agriculture,  hydrology,  and  meteorology.  In  the  field 

This  report  follows  the  style  and  format  of  the  Journal  of  Climate  and 
Appl led  Meteorology. 
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of  agriculture  Idso  et  al_.  (1975a)  described  the  need  for  soil  moisture 
information  for  many  diverse  applications,  including  improved  crop 
yield  forecasting  and  irrigation  scheduling.  In  fact,  knowledge  of  the 
moisture  condition  of  the  soil  is  needed  before  undertaking  any  agri¬ 
cultural  activity  from  planting  to  harvesting. 

Remote  sensing  has  been  used  to  evaluate  or  monitor  numerous 
hydrologic  parameters.  Burgy  (1974)  cites  thirteen  hydrologic  parame¬ 
ters  which  can  be  measured  or  inferred  through  remote  sensing.  These 
are  surface  temperature,  cloud  cover,  atmospheric  liquid  water  content, 
snow  cover,  water  equivalent  of  snowpack,  vegetation,  land  use,  soil 
moisture,  short  wave  radiation  in  the  radiation  budget,  long  wave 
radiation  in  the  radiation  budget,  water  body  configuration,  soil 
types,  and  basin  area.  Knowledge  of  the  moisture  content  of  the  soil 
surface  layer  is  necessary  for  separating  rainfall  into  its  runoff  and 
infiltration  components.  This  knowledge  is  needed  for  both  soil  mois¬ 
ture  modeling  and  flood  forecasting. 

In  meteorology,  soil  moisture  determines  the  division  of  net  radi¬ 
ation  into  its  latent  and  sensible  heat  components.  An  estimation  of 
evaporation  or  evapotranspiration  requires  input  of  current  soil  mois¬ 
ture  conditions.  Most  radiation,  water  budget,  or  global  circulation 
models  need  an  estimate  of  these  parameters  as  initial  input. 

Passive  Microwave  Indication  of  Soil  Moisture 

Emitted  radiation  at  microwave  frequencies  is  a  function  of  the 
moisture  content  of  the  emitting  soil  layer.  This  relationship  has 


been  demonstrated  by  Investigations  by  Cihlar  and  Ulaby  (1975), 
Schmugge  (1977),  Schmugge  et_  £h  (1974,  1976,  and  1978),  and  Newton 
(1977).  Basically,  this  relationship  is  a  function  of  the  extremely 
high  dielectric  constant  of  water,  and  the  relatively  low  values  for 
air  and  dry  soil.  Since  the  dielectric  constant  inversely  determines 
the  soil  emissivity,  an  increase  in  volumetric  water  content  causes  a 
decrease  in  the  emissivity.  The  brightness  temperature,  or  radiation 
received  at  the  sensor  antenna,  is  equal  to  the  emissivity  times  the 
temperature  of  the  emitting  layer.  It  is  possible  to  estimate  the 
emissivity  and  therefore  the  soil  moisture  content  by  measuring  the 
brightness  temperature  and  estimating  the  soil  surface  temperature  from 
meteorological  observations.  Meneely  (1977)  attempted  to  show  this  by 
using  1.55  cm  passive  microwave  data  from  the  Nimbus  5  sensor,  however 
he  did  not  get  particularly  encouraging  results.  Later  work  by 
McFarland  and  Blanchard  (1977),  fheis  (1979),  and  Blanchard  et  al . 
(1980,  1981)  involving  1.55  cm  passive  microwave  data.  Indicated  a 
strong  linear  correlation  with  rainfall  history  and  inferred  soil  mois¬ 
ture  for  essentially  bare  agricultural  soils  through  a  time  series 
analysis  of  brightness  temperatures.  Harder  (1984)  found  similar  cor¬ 
relations  using  Nimbus  7  passive  microwave  data  at  the  1.66  cm  wave¬ 
length. 


Ground  Observations 

One  of  the  largest  problems  confronting  all  researchers  is  obtain¬ 
ing  accurate  and  reliable  ground  observations  (Harder,  1984).  In 
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satellite  remote  sensing  the  vast  area  sensed  presents  an  even  greater 
problem  in  acquiring  accurate  ground  observations.  There  are  many  con¬ 
ventional  ways  to  measure  soil  moisture  over  limited  areas,  but  a  prac¬ 
tical  method  for  measuring  soil  moisture  over  an  area  consisting  of 
several  thousands  of  square  kilometers  has  not  been  developed.  There¬ 
fore  this  investigation  will  use  the  Antecedent  Precipitation  Index 
(API)  to  infer  the  ground  observations  needed  to  estimate  soil  mois¬ 
ture.  This  index  represents  soil  moisture  as  a  first  order  Markov 
chain  through  input  of  effective  precipitation  (Harder,  1984). 


Analytical  Approach 


This  investigation  involved  both  a  qualitative  and  quantitative 
approach.  Initially,  the  data  were  analyzed  qualitatively  to  determine 
all  patterns  or  trends.  The  qualitative  analysis  consisted  of  a  tem¬ 
poral  analysis  of  SMMR  brightness  temperature  data  between  25  Oct  1978 
and  10  Nov  1979,  areal  analysis  on  one  case  study  grid  cell,  and  image 
processor  files  for  all  grid  cells  in  the  study  area.  One  day  was 
chosen  from  the  time  series  to  be  used  as  a  case  study.  Areal  maps  of 
passive  microwave  brightness  temperatures  over  the  study  area  were  used 
as  input  to  create  various  image  processor  files.  The  visual  patterns 
evident  on  the  image  processor  files  were  then  simultaneously  compared 
to  determine  which  microwave  channel  had  the  best  qualitative  relation¬ 
ship  to  surface  soil  moisture  conditions.  By  knowing  the  actual  sur¬ 
face  moisture  condition  through  supplementary  data,  one  can  determine 
qualitatively  which  channel  of  microwave  data  has  the  best  visual 
relationship  to  the  surface  moisture  condition. 
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The  quantitative  approach  involved  regression  analysis  relating 
the  passive  microwave  data,  along  with  various  combinations  of  the  pas¬ 
sive  microwave  brightness  temperatures,  with  the  API.  The  microwave 
data  was  analyzed  separately  in  terms  of  polarization  and  monthly, 
bi-monthly,  and  annual  response.  This  analysis  will  determine  which 
wavelength  and  polarization  combination  gives  the  best  indication  of 
the  soil  moisture,  as  inferred  by  the  API. 
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CHAPTER  II 

LITERATURE  REVIEW 

Passive  Microwave  Remote  Sensing  of  Soil  Moisture 

There  are  numerous  methods  using  microwave  measurements  to  esti¬ 
mate  soil  moisture.  These  methods  are  sunmarized  by  Idso  et  ^1_. 
(1975a,  1975b),  Jackson  et_  aL  (1978),  Schmugge  et  (1976,  1980), 
Schmugge  (1978,  1980a,  1980b),  Schmugge  and  Choudhury  (1981),  Newton 
(1980),  and  Theis  et^  aU  (1982).  The  microwave  methods,  both  active 
and  passive,  deal  with  the  1-50  cm  wavelength  region.  These  methods 
rely  on  the  dielectric  properties  of  the  soil  and  water  mixture.  Since 
the  dielectric  properties  of  the  substance  determine  the  electromag¬ 
netic  wave  propagation  characteristics  within  the  substance,  these 
properties  will  also  affect  the  surface  emissive  and  reflective  proper¬ 
ties.  The  microwave  emission  and  reflection  are  then  a  direct  function 
of  the  soil  moisture  and  can  be  measured  by  active  (radar)  scattero- 
meters  or  passive  radiometers.  Active  microwave  methods  are  discussed 
extensively  in  Ulaby  ^t  a]_.  (1974),  Ulaby  (1974),  Ulaby  and  Batlivala 
(1974),  Batlivala  and  Ulaby  (1977),  Ulaby  et_  al_.  (1977,  1978)  and 
Schmugge  et  ^1_.  (1979). 

Passive  microwave  radiometers  measure  the  surface  thermal  emis¬ 
sion.  In  the  microwave  region,  the  Raylelgh-Jeans'  approximation  shows 
that  the  Intensity  of  the  observed  emission,  called  brightness  tempera¬ 
ture,  is  essentially  proportional  to  the  product  of  the  surface  temper¬ 
ature  and  the  emissivity.  Advantages  of  passive  microwave  approaches 
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Include  the  ability  to  penetrate  non-precipitating  clouds  and  light 
vegetation.  A  few  disadvantages  are  the  sensor  response  only  involves 
the  top  0-2  cm  of  the  soil,  and  relatively  coarse  data  resolution  (>10 

km),  as  a  function  of  the  altitude  and  diameter  of  the  antenna  of  the 

radiometer.  The  formula  for  the  brightness  temperature  involves  the 
physics  of  passive  microwave  remote  sensing.  Since  the  main  emphasis 
of  this  investigation  will  involve  the  passive  microwave  region,  the 
remote  sensing  characteristics  of  this  region  will  be  discussed  more 

throughly  in  the  next  section. 

Formula  for  Passive  Microwave  Emission 

According  to  Lintz  and  Simonett  (1976),  microwave  emission  from  an 
object  can  be  described  best  by  the  concept  of  a  theoretical  blackbody, 
a  non-reflecting  opaque  object.  Planck's  law  gives  an  expression  for 
the  spectral  radiant  emittance  as  a  function  of  the  wavelength  for 

black-body  radiation  at  a  given  temperature.  The  formula  is 

>fl  <» 

where 

Mx  =  spectral  radiance  of  power  emitted  (watts  cm" 2ym"  l) 
h  *  Plank's  constant  (6.626xl0"34  J  sec) 
k  *  Boltzman's  constant  (1.380xl0-23  JK"1) 

A  3  wavelength  (m) 

T  *  blackbody  temperature  (K) 
c  »  velocity  of  light  (2.998x10s  m  sec"1) 
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Planck's  law  for  an  ideal  blackbody  in  the  microwave  region  is 
given  by  the  Rayleigh-Jeans’  approximation  of  Planck's  radiation 
formula.  This  approximation  can  be  derived  by  first  taking  a  Maclaurin 
series  expansion  of  the  bracketed  quantity  in  (1).  A  general  form  for 
this  expansion  is 

eX  =  1  +  x  +  7T  +  TT  +  •  (2 ) 


Expanding  the  exponent  in  this  manner  yields 


he  2/  xkT  ,  he  3/  xkT  . 

~v. — +  — + 


(3) 


At  microwave  frequencies  he  «  kT,  therefore,  we  can  neglect  the 
second  and  all  successive  terms  in  (3).  This  equation  then  reduces  to 
the  Rayleigh-Jeans'  approximation. 


M 


X 


(4) 


This  gives  the  total  spectral  radiance.  The  polarized  radiance, 
where  each  plane  of  polarization  contains  one-half  of  the  total 
radiance,  would  be  given  by 


M 


P 


—r  =  e  oT"  (for  non-microwave  regions) 
x^ 


(5) 


where  e  is  the  emissivity  and  o  is  the  Stefan-Boltzmnn  constant, 
5.67xl0-12  w  cm-2K-4. 
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At  microwave  wavelengths,  the  spectral  radiance  is  linearly  pro¬ 
portional  to  the  temperature  of  the  emitting  layer.  At  shorter  wave¬ 
lengths,  the  spectral  radiance  is  proportional  to  the  fourth  power  of 
the  temperature.  This  linear  relationship  has  led  to  the  expression  of 
the  spectral  radiance  as  a  brightness  temperature.  For  non-blackbod- 
ies,  the  brightness  temperature  is  equal  to  the  product  of  the  emis- 
sivity  and  the  temperature  of  the  emitting  layer.  The  radiance  of  an 
equivalent  blackbody  is  multiplied  by  a  relative  efficiency  factor 
called  the  emissivity.  Equation  (5)  can  be  rearranged  to  get  the 
polarized  brightness  temperature  (Tb),  defined  by 


which  reduces  to 

Tb  =  eT  (7) 


where  e  is  the  emissivity  and  T  is  the  graybody  temperature. 

The  Tb  measured  by  a  satellite  radiometer  is  given  by  Schmugge 
(1980a)  as 


Tb  =  Tr[rTsky  +  (l-r)Tsoil]  +  Tatm  (8) 

where  Tr  is  the  atmospheric  transmission  and  r  is  the  surface 
reflectivity.  This  formula  separates  the  Tb  into  its  three  compo¬ 
nents.  These  are  the  reflected  sky  brightness  temperature  (Tsky),  the 
thermal  soil  emission  or,  under  conditions  of  non-transmittance  and 


thermal  equilibrium,  emissivity  (c,  since  e  =  1-r),  and  the  contribu¬ 
tion  due  to  the  intervening  atmosphere  (Tatm).  The  reflected  sky 
brightness  temperature  is  the  product  of  the  downwelling  brightness 
temperature  and  the  surface  reflectivity,  minus  the  loss  due  to  the 
intervening  atmosphere.  The  sky  radiation  or  sky  brightness  tempera¬ 
ture  consists  of  the  4K  celestial  background  emission  plus  downwelling 
atmospheric  radiation. 

Weger  (1960)  computed  the  sky  brightness  temperature  for  numerous 
sky  conditions  for  wavelengths  between  0.43  and  3.00  cm.  Sky  bright¬ 
ness  values  for  a  representative  wavelength  of  1.80  cm  at  a  viewing 
angle  of  50  degrees,  were  about  10  K  for  clear  skies,  15  K  for  moderate 
cloud,  and  30  K  for  uniform  moderate  precipitation.  Under  clear  skies, 
the  atmospheric  transmission  will  be  close  to  1.0.  With  a  high  surface 
reflectivity  (90%),  the  total  contribution  will  be  close  to  9  K.  Under 
precipitating  conditions,  the  atmospheric  transmission  and  surface 

reflectivity  will  decrease.  With  an  atmospheric  transmission  of  80% 
and  a  surface  reflectivity  of  30%,  the  contribution  will  be  7.2  K. 

This  shows  that  the  maximim  component  of  the  Tb  due  to  the  reflected 
sky  brightness  temperature  could  be  as  large  as  9  K. 

The  thermal  emission  from  the  soil  is  generated  within  the  soil 

layer.  The  amount  of  energy  generated  at  a  particular  point  depends  on 
the  soil  temperature  at  that  point  and  the  volumetric  dielectric 

properties.  As  the  energy  propagates  from  the  point  of  origin,  it  is 
modified  by  the  dielectric  gradient  it  encounters.  Numerous  radiative 
transfer  models  have  been  developed  to  account  for  these  energy 
exchanges.  Some  of  these  models  are  described  in  Wang  and  Schmugge 
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(1978,  1980),  Burke  _eta]_.  (1979),  and  Schmugge  and  Choudhury  (1981). 
Radiative  Transfer  Models 

Most  of  the  radiative  transfer  models  require  a  fairly  detailed 
knowledge  about  the  soil  mediun  and  some  knowledge  about  the  surface 
parameters.  One  very  significant  result  from  these  models  was  given  by 
Wilheit  (1978).  Wilheit  defined  thermal  and  reflective  sampling  depths 
and  showed  that  the  thermal  depth  is  on  the  order  of  one  wavelength, 
while  the  reflective  depth  was  on  the  order  of  a  few  tenths  of  a  wave¬ 
length.  The  thermal  sampling  depth  referred  to  the  layer  whose  effec¬ 
tive  temperature  determined  the  amount  of  energy  available  for  emis¬ 
sion.  The  reflective  sampling  depth  referred  to  the  layer  whose  mois¬ 
ture  content  determined  the  emissivity  of  the  soil.  Wang  et^  aK 
(1980a)  used  the  empirical  model  of  Wang  and  Schmugge  (1980)  to  deter¬ 
mine  the  emission  from  the  soil  surface.  They  computed  the  expected 
normalized  brightness  temperature  or  effective  emissivity,  for  both  6 
and  21  cm  wavelengths,  of  a  soil  surface  covered  by  different  types  of 
vegetation.  While  both  observations  and  predictions  showed  a  general 
decrease  of  sensitivity  with  increasing  vegetation,  as  theory  suggests, 
they  did  notice  some  discrepancies  between  measurements  and  theory. 
They  suggested  that  the  discrepancies  may  be  explained  by  roughness, 
which  was  not  specifically  accounted  for  in  the  model. 

The  effects  due  to  the  Intervening  atmosphere  fall  into  two 
categories,  the  direct  radiation  from  the  atmosphere  incident  upon  the 
sensor,  and  the  atmospheric  attenuation  of  reflected  or  emitted  surface 
radiation.  At  long  wavelengths  these  effects  are  minimal  when  compared 


to  the  effects  from  the  other  components  and  can  be  neglected  in  a 
first  approximation  of  soil  moisture. 
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Simplified  Soil  Emission  Models 

The  radiative  transfer  models  described  in  the  literature  often 
require  extensive  knowledge  of  the  soil  structure,  surface  roughness, 
vegetation,  and  a  wide  continuun  of  meteorological  elements.  These 
detailed  data  are  extremely  difficult  to  obtain  over  the  large  areas 
normally  covered  by  low-resolution  satellites.  Most  simplified  models 
represent  the  soil  as  a  single  emitting  layer.  This  makes  the  models 
well  suited  for  using  data  from  satellite-based  passive  microwave 
radiometers  whose  depth  of  penetration  is  only  a  few  centimeters. 

.•se  models  rely  on  the  strong  dependence  of  the  soil  dielectric  prop¬ 
erties  on  the  soil  moisture  content.  The  magnitude  of  this  dependence 
is  so  great  that  it  overshadows  all  other  components  of  the  complex 
radiative  transfer  models.  These  secondary  components  can  therefore  be 
estimated  or  even  neglected  making  a  simplified  model  applicable  to 
satellite  remote  sensing.  One  of  the  most  important  results  from  one 
of  these  simplified  models  was  given  by  Newton  (1977).  Newton  present¬ 
ed  model  results  showing  that  brightness  temperature  was  a  nearly 
linear  function  of  soil  moisture  content  throughout  the  moisture 
content  range,  except  near  the  dry  end.  The  lack  of  a  linear  relation¬ 
ship  at  the  dry  end  was  explained  by  the  tight  molecular  bond  between 
the  water  and  soil  molecules  in  this  moisture  range.  He  also  observed 
that  when  the  viewing  angle  was  increased  from  nadir  to  40  degrees,  the 
vertically  polarized  brightness  temperature  component  increased  and  the 


horizontally  polarized  brightness  temperature  component  decreased.  He 
stated  that  the  mean  of  the  two  polarized  components  was  nearly  inde¬ 
pendent  of  viewing  angle.  This  is  shown  by  Fig.  1. 


Wang  and  Schmugge  (1978,  1980)  used  a  simple  empirical  model  of 
the  dielectric  behavior  of  soil-water  mixtures.  This  was  a  two-layer 
model  with  an  adjustable  lower  layer.  They  determined  that  the  dielec¬ 
tric  constant  of  a  soil  was  not  very  responsive  to  soil  moisture  below 
a  certain  transition  moisture  value,  and  that  this  transition  value  was 
a  simple  linear  function  of  the  permanent  wilting  point  of  the  soil. 
Since  the  transition  moisture  value  was  dependent  on  soil  texture,  soil 
texture  maps  would  be  useful  in  predicting  the  soil  dielectric  behav¬ 
ior.  Wang  (1980)  expanded  on  this  observation  when  he  stated  that 
wavelengths  shorter  than  30  cm  have  a  moisture  response  that  depends  on 
soil  texture,  while  longer  wavelengths  have  no  such  dependence. 

Noise 

The  relationship  between  the  soil  moisture  signal  and  the  sensor 
response  is  affected  by  various  noise  factors.  These  noise  factors 
include  vegetation,  roughness,  and  instrument  errors.  The  most  signif¬ 
icant  noise  factors  in  the  passive  microwave  region  are  vegetation  and 
roughness.  The  passive  microwave  response  will  be  a  composite  of  the 
emission  from  the  total  area  sensed.  This  response  even  over  "bare" 
agricultural  soils  will  contain  vegetative  effects  resulting  from 
trees,  pastures,  or  roadside  ditches.  The  manner  in  which  roughness 
and  vegetation  affect  the  soil  moisture  signal  is  different  but  the 
result  is  the  same.  Both  factors  decrease  the  soil  moisture  signal 
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Figure  1 
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Angular  dependence  of  the  polarized  brightness 
temperature.  A  schematic  representation  of  the 
effect  of  soil  moisture  on  brightness  temperature 
and  polarization  {adapted  from  Harder.  1984). 


detectable  at  the  sensor.  Poe  and  Edgerton  (1971)  investigated  respon¬ 
ses  from  airborne  microwave  radiometers  at  wavelengths  of  0.32,  0.81, 
1.55,  6.0,  and  21  cm.  They  concluded  that  the  shortest  three  wave¬ 
lengths  were  too  strongly  influenced  by  noise,  particularly  roughness, 
to  be  of  much  use  in  soil  moisture  estimation.  The  longest  two  wave¬ 
lengths  were  consistently  sensitive  to  soil  moisture  even  in  the 
presence  of  a  vegetative  canopy.  Choudhury  et  _al_.  (1979)  presented 
model  calculations  which  agreed  with  previous  experimental  results,  and 
concluded  roughness  effects  were  greatest  for  wet  soils.  Wang  et  al . 
(1980b)  found  that  vegetation  degrades  soil  moisture  sensitivity  at 
both  6  and  21  cm  wavelengths  by  an  amount  related  to  the  amount  of 
biomass  and/or  moisture  content  of  the  vegetation.  Wang  and  Choudhury 
(1981)  developed  a  microwave  emission  model  that  had  a  two-parameter 
representation  of  soil  surface  roughness.  They  suggested  both  the 
polarization  mixing  parameter  and  the  roughness  height  parameter  could 
be  estimated  empirically  by  a  time  series  analysis  of  transformations 
of  polarized  passive  microwave  brightness  temperatures.  McFarland 
(1982)  suggested  that  when  vegetation  obscures  the  soil  moisture 
signal,  the  observed  brightness  temperature  would  indicate  crop  temper¬ 
ature.  This  crop  temperature  could  be  compared  to  the  ambient  tempera¬ 
ture  to  give  an  indication  of  crop  stress.  Inferred  values  of  crop 
stress  could  then  be  used  to  infer  qualitatively  soil  moisture  content 
in  the  root  zone. 

Instrument  errors  including  calibration  cause  an  additional  source 
of  noise  in  the  passive  microwave  brightness  temperature  measurements. 
All  Instruments  have  an  associated  error  with  their  measurements. 
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however  this  error  is  greatly  reduced  by  routine  calibration.  If  the 
sensor  is  not  calibrated  and  checked  for  consistency,  the  data  becomes 
suspect.  While  there  is  an  error  range  associated  with  the  measure¬ 
ments  from  these  instruments,  it  is  very  small  when  compared  to  the 
other  sources  of  error  such  as  vegetation  and  roughness,  and  the  inter¬ 
vening  atmosphere. 

McFarland  and  Blanchard  (1977)  and  B1  anchard  et  al_.  (1981)  report¬ 
ed  data  interpretation  problems  due  to  the  cross-track  scan  pattern  of 
the  Nimbus  5-ESMR  sensor  which  resulted  in  a  non-uniform  viewing  angle 
and  a  degrading  resolution  cell  as  the  viewing  angle  increased.  Even 
though  the  Nimbus  7  SMMR  has  a  constant  viewing  angle,  the  data  still 
suffer  from  a  degrading  resolution  cell,  like  the  Nimbus  5-ESMR  data, 
as  the  wavelength  increases.  While  the  magnitude  of  the  instrument 
errors  from  the  Nimbus  5-ESMR  was  significant,  effects  of  vegetation 
and  roughness  will  totally  mask  any  discernible  errors  from  the  Nimbus 
7  SMMR  sensor. 


Soil  Moisture  Models 

Numerous  soil  moisture  models  are  available  in  the  literature 
ranging  from  extremely  complex  models  requiring  considerable  meteor¬ 
ological  input  along  with  detailed  soil  and  vegetation  parameters,  to 
relatively  simple  models  which  can  utilize  low-resolution  satellite 
data. 

Estimation  of  Evapotranspiration 


One  of  the  largest  problems  encountered  in  developing  a  soil 
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moisture  model  is  the  method  used  to  estimate  evaporation  or 
evapotranspi ration.  Penman  (1956)  gave  some  interesting  quotes  on  this 
subject  from  previous  authors.  Symons  laid  the  groundwork  for  future 
studies  in  this  area  when  he  said  in  1867  “evaporation  is  the  most 
desperate  branch  of  the  desperate  science  of  meteorology 
Maluschitski  followed  this  example  by  stating  in  1900  that  "no 
correlation  can  be  established  between  the  evaporation  from  a  water 
surface  and  that  for  a  cultivated  soil,  and  still  less  in  the  case  of  a 
soil  covered  with  plants."  Cleveland  Abbe  continued  this  trend  in  1905 
when  he  said  “of  course  (we)  need  to  know  the  loss  of  water  by 
evaporation  but  in  nature  this  is  so  much  mixed  up  with  seepage, 
leakage  and  consumption  by  ...  plants  that  our  meteorological  data  are 
of  comparatively  little  importance."  We  have  progressed  forward  since 
these  statements  were  made,  however  a  unique  method  for  estimating 
evaporation  or  evapotranspiratlon  still  has  not  been  developed. 

According  to  Schmugge  (1978)  some  of  the  major  differences  in  the 
models  involve  the  method  of  computation  or  estimation  of  potential 
evaporation,  how  infiltration  or  runoff  is  computed,  the  way  precipita¬ 
tion  or  evaporative  demand  is  defined,  the  method  of  computing  soil 
evaporation  and  plant  transpiration,  the  number  of  soil  layers  to  use 
and  whether  they  are  saturated  or  unsaturated  layers,  and  the  way  the 
thermal  properties  of  the  soil  system  are  accounted  for.  The  major 
advantage  of  using  a  model  is  to  provide  timely  point  or  areal  esti¬ 
mates  of  soil  moisture  over  areas  where  these  estimates  are  not  readily 
available.  The  major  disadvantage  is  the  inherent  Inaccuracy  of  the 
estimate. 
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Chang  (1968)  classified  the  various  methods  to  estimate  evapora¬ 
tion  or  evapotranspi ration  into  five  major  classes.  These  are  direct 
measurements,  aerodynamic  methods,  energy  balance  approaches,  empirical 
methods,  and  evaporation  pans  or  atmometers.  The  input  to  soil  mois¬ 
ture  models  normally  involves  one  of  the  nunerous  empirical  methods  by 
using  an  empirical  moisture  depletion  coefficient.  If  this  coefficient 
is  fixed,  then  the  model  is  applicable  only  to  areas  with  similar 
environments  to  the  one  in  which  it  was  developed.  Often  the  model 
will  have  only  limited  qualitative  usefulness  due  to  inherent  errors, 
even  over  similar  regions.  There  seems  to  be  a  model  in  the  literature 
for  nearly  every  possible  situation,  therefore  only  a  few  representa¬ 
tive  soil  moisture  models  will  be  discussed  here. 

API  Model 

The  API  is  a  very  simple  soil  moisture  model  that  can  rely  predom¬ 
inantly  on  passive  microwave  data.  The  resultant  soil  moisture  values 
are  much  coarser  than  detailed  models  due  to  the  low- resolution  input 
data  and  the  greatly  simplified  physics  of  the  model.  Since  the  prob¬ 
lem  of  obtaining  an  accurate  and  unique  model  for  estimating  evapora¬ 
tion  or  evapotranspi  ration  has  not  been  solved,  this  method  may  prove 
to  be  a  viable  alternative. 

Saxton  and  Lenz  (1967)  used  a  modified  API  which  they  termed  an 
antecedent  retention  index.  They  used  an  estimate  of  the  potential 
evaporation  to  deplete  moisture  during  the  constant  rate  stage  of 
drying,  and  an  exponential  depletion  rate  in  the  falling  rate  stage. 
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An  excellent  description  of  these  drying  stages  is  given  by  Derendinger 
(1971). 

According  to  Li nsley  _et  _aK  (1975)  the  most  common  soil  moisture 
index  used  for  storm  runnoff  models  is  the  API  in  the  form  used  by 
Kohler  and  Li nsley  (1951).  This  form  of  the  API  was 

It  =  IoK  (9) 

where  Io  is  the  initial  value,  K  is  the  depletion  coefficient, 
less  than  one,  and  t  is  the  time  in  days.  They  stated  that  the  deple¬ 
tion  coefficient  or  recession  factor  should  be  a  function  of  potential 
evaporation  or  season.  They  suggested  a  seasonal  value  could  be  used 
since  it  would  be  strongly  related  to  potential  evapotranspi ration. 
They  noted  also  that  initialization  errors  disappear  after  a  few  weeks 
since  the  errors  decayed  logarithmically. 

Choudhury  and  Blanchard  (1981)  used  a  water  balance  approach 
taken  from  Global  Climate  Models  to  derive  several  forms  of  the  API  as 
a  first-order  Markov  chain.  Their  formula  for  the  recession  factor  was 

k  =  exp  t-E/Z(FC-WP)]  (10) 

where  E  is  the  daily  total  potential  evapotranspi  ration  (PET)  for  the 
previous  day,  Z  is  the  soil  thickness,  FC  is  the  field  capacity  of  the 
soil,  and  WP  is  the  permanent  wilting  point.  They  noted  that  a  first 
approximation  of  the  recession  factor  could  be  obtained  from  soil  tex¬ 
ture  information  and  a  climatological  estimate  of  PET. 

The  form  of  the  API  used  by  McFarland  and  Harder  (1982)  and 
Harder  (1984)  for  use  with  passive  microwave  studies  was 
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API  (i )  =  API (1-1)  *  k  +  P(i )  (11) 

where  API ( i )  is  the  API  value  for  day  i,  API(i-l)  is  the  previous  day's 
value,  k  is  a  depletion  coefficient  less  than  one,  and  P(i)  is  the 
effective  precipitation  for  day  i.  The  months  of  February  and  August 
were  chosen  to  make  the  depletion  coefficient  a  simple  sinusoidal  model 
which  approximated  the  annual  temperature  minimum  and  maximim  values. 
McFarland  (1976)  initially  used  satellite  measurements  to  infer  soil 
moisture  through  correlations  between  these  measurements  and  API.  He 
used  a  constant  k  value  of  0.9  for  the  June  and  August  correlations 
with  the  L-band  (21  cm)  brightness  temperatures  from  the  Skylab  mis¬ 
sion.  McFarland  and  Blanchard  (1977)  used  a  varying  k  factor  as  a 
function  of  season  in  order  to  remove  the  seasonal  temperature  trend 
effects  from  the  three  month  period  of  record.  B1  anchard  et  al_.  (1980) 
stated  that  the  general  form  of  a  curve  which  represented  inverted  mean 
daily  temperature,  represented  seasonal  changes  in  k.  This  curve  of 
inverted  mean  daily  temperature  was  found  to  be  remarkably  similar  to 
the  inverse  of  a  curve  of  the  mean  monthly  potential  evapotranspi ra¬ 
ti  on,  calculated  by  the  Penman  equation.  Making  the  depletion  coef¬ 
ficient  a  function  of  season  agreed  with  previous  work  by  Kohler  and 
Linsley  (1951).  The  effective  precipitation  was  computed  by  raising 
precipitation  to  a  power  of  0.891  as  suggested  by  Blanchard  et  ^1_. 
(1980).  McFarland  and  Harder  (1982)  applied  this  model  by  using  simple 
linear  regression  to  correlate  the  API  with  passive  microwave 
emissivity  values.  The  emissivities  were  computed  from  the  ratio  of 
the  microwave  brightness  temperatures  to  the  air  temperature. 
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Other  Models 


Baier  and  Robertson  (1966)  developed  a  multi-layered  Versatile 
Budget  (VB)  model  similar  to  the  API  model,  but  with  a  more  detailed 
depletion  coefficient.  This  coefficient  was  computed  from  the  product 
of  the  potential  evapotranspiration,  an  adjustment  for  a  specific  dry¬ 
ing  curve,  and  a  soil-plant  characteristic  coefficient.  The  previous 
day's  soil  moisture  value  was  then  multiplied  by  the  depletion  coef¬ 
ficient  to  get  the  effective  precipitation  to  be  added  to  the  infiltra¬ 
tion  amount. 

Van  Bavel  and  Lascano  (1980)  developed  a  model  called  CONSERVB 
that  calculated  soil  moisture  and  temperature  profiles  from  weather 
observations.  This  energy  balance  approach  required  either  hourly  or 
daily  observations  of  global  radation,  wind  speed,  air  temperature, 
humidity,  and  precipitation  amount  and  duration.  These  are  the  same 
parameters  required  for  the  widely-used  Penman  potential  evapotrans¬ 
piration  equation.  Smith  and  Newton  (1983)  modified  the  CONSERVB  model 
to  accept  microwave  radiometer  data  instead  of  meteorological  observa¬ 
tions.  They  used  data  simulation  for  both  6  and  21  cm  wavelengths  to 
compare  to  the  original  CONSERVB  model.  The  resulting  soil  moisture 
predictions  were  nearly  identical  to  the  original  model.  This  experi¬ 
ment  shows  that  microwave  data  can  be  substituted  for  actual  meteoro¬ 
logical  observations  with  an  insignificant  loss  of  predictive  capabil- 


Surface  and  Aircraft  Microwave  Investigations 


Surface  Emission 

Edgerton  _et  _a1_.  (1971)  measured  the  passive  microwave  emission 
from  several  substances  including  soils,  sediments,  sea  ice,  snow,  and 
snowpacks.  Combinations  of  viewing  angle  and  wavelength  were  used  to 
examine  different  types  of  soil  texture.  They  concluded  that  the 
brightness  temperature  was  a  function  of  radiometer  polarization  and 
wavelength,  thermodynamic  temperature,  soil  density,  surface  roughness, 
vegetative  cover,  and  soil  moisture  content.  The  variation  of  bright¬ 
ness  temperature  with  soil  texture  and  viewing  angle  was  inherently 
accounted  for  in  the  experimental  design.  They  found  that  the  2.2  cm 
wavelength  response  was  related  strongly  to  both  soil  moisture  and  the 
soil  temperature  of  the  top  centimeter. 

Kirdyashev  et  aU  (1979)  undertook  a  similar  study  when  they 
investigated  the  sensitivity  of  the  microwave  emission  under  varying 
types  and  amounts  of  vegetation.  They  analyzed  the  radiometric 
response  at  several  viewing  angles  at  wavelengths  ranging  from  0.8  to 
30  cm.  They  found  that  the  sensitivity  of  the  microwave  emission 
decreased  as  the  amount  of  vegetation  increased,  and  presented  quanti¬ 
tative  models  of  the  vegetation  attenuation  effects.  They  suggested 
the  possibility  of  using  microwave  data  to  estimate  biomass,  plant 
height,  and  water  content. 

Investigations  Involving  Polarized  Data 


Newton  (1977)  showed  that  increasing  the  sensor's  viewing  angle 
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from  nadir  to  near  50  degrees  caused  an  increase  in  the  vertically 
polarized  brightness  temperature  component  and  a  decrease  in  the 
horizontally  polarized  brightness  temperature  component.  He  showed 
that  the  mean  of  the  two  polarized  components  was  nearly  independent  of 
viewing  angle.  When  discussing  the  effects  of  vegetation  and  surface 
roughness,  Newton  showed  that  the  effect  of  both  of  these  parameters 
was  to  decrease  the  sensitivity  of  the  emitted  brightness  temperature 
to  soil  moisture.  He  concluded  that  the  difference  between  the  L-band 
vertically  polarized  component  of  the  brightness  temperature  and  the 
L-band  horizontally  polarized  component  at  a  viewing  angle  of  35 
degrees  could  be  used  to  estimate  surface  roughness. 

Burke  and  Paris  (1975)  and  Burke  et_  ^1_.  (1979)  also  examined  a 

combination  of  the  polarized  components  of  the  brightness  tempera¬ 
tures.  They  investigated  the  sensitivity  of  the  dual  polarized  micro- 
wave  response  to  soil  moisture  at  wavelengths  of  2.8  and  21  cm.  They 
compared  values  of  the  first  two  Stokes'  parameters,  P  =  0.5(Tv+Th), 
and  Q  =  (Tv-Th),  and  found  they  could  identify  both  soil  moisture  and 
surface  roughness  effects.  Q  was  found  to  be  very  sensitive  to  soil 
moisture  at  gravimetric  moisture  contents  less  than  15  percent,  while  P 
was  more  sensitive  at  moisture  contents  greater  than  15  percent.  Q  was 
more  closely  related  with  surface  roughness  and  nearly  invariant  with 
temperature,  while  P  was  found  to  be  nearly  invariant  with  viewing 
angle. 

Burke  (1980)  continued  the  evaluation  of  the  Stokes'  parameters 
using  aircraft  data  at  the  Ku,  X,  and  L-band  wavelengths.  She  found 
that  the  first  Stokes'  parameter,  P  =  0.5(Tv+Th),  was  more  sensitive  to 
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soil  moisture  than  either  of  the  individual  components  of  the  polarized 
brightness  temperature.  The  second  Stokes'  parameter,  Q  =  (Tv-Th), 
showed  considerable  scatter  when  related  to  surface  roughness.  This 
disagreed  with  her  previous  work,  where  Q  had  been  found  to  be  more 
closely  related  with  surface  roughness  than  P.  She  agreed  with  pre¬ 
vious  observations  that  vegetation  completely  obscurred  any  response  to 
soil  moisture  from  the  shorter  wavelengths,  and  that  the  L-band  wave¬ 
length  was  relatively  insensitive  to  vegetation.  She  concluded  that 
multi-spectral,  multi-temporal  microwave  observations  could  indicate 
shallow  subsurface  moisture  profiles. 

Wang  _et  ^1_.  (1981)  took  the  ratio  of  the  two  Stokes'  parameters 

and  defined  it  as  a  polarization  factor  (P),  where  P  =  [Tv-Th/0.5 
(Tv+Th)].  They  used  this  factor  to  evaluate  the  results  of  ground- 
based  measurements  of  brightness  temperatures  at  wavelengths  of  6 
(C-band)  and  21  (L-band)  cm.  They  found  that  the  ratio  of  C-band  to 
L-band  brightness  temperatures  could  discriminate  between  bare  and 
vegetated  fields.  In  addition  they  found  that  P  could  indicate  bio¬ 
mass,  with  high  values  of  P  corresponding  to  bare  fields  or  dry  corn 
fields. 

Normalized  Brightness  Temperatures 

Mo  .££.  (1980)  and  Mo  and  Choudhury  (1980)  showed  that  soil 

moisture  was  related  to  an  emissivity  or  a  normalized  brightness  tem¬ 
perature  (Tb),  where  Tb  was  defined  as  brightness  temperature  divided 
by  the  effective  temperature  of  the  emitting  layer.  They  showed  that 
the  effective  temperature  was  approximately  equal  to  the  surface 
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temperature  for  short -wavelength  sensors,  but  that  the  entire  thermal 
sampling  depth  must  be  considered  for  long-wavelength  sensors.  Wang  et 
a] .  (1981)  used  this  definition  for  normalized  brightness  temperatures 

when  they  analyzed  microwave  data  obtained  at  wavelengths  of  2.8,  6.0, 
21,  and  50  cm.  They  defined  the  emitting  layer  as  the  soil  layer, 
whose  moisture  content  was  being  measured,  with  a  sampling  depth  equal 
to  about  one  tenth  of  a  wavelength.  They  found  that  the  wavelength 
with  the  greatest  moisture  sensitivity  was  6.0  cm  for  bare  soil,  21  cm 
for  alfalfa,  and  50  cm  for  grass. 

Njoku  and  O'Neill  (1982)  also  used  a  normalized  brightness  temper¬ 
ature  in  an  investigation  of  the  relationship  of  the  volumetric  soil 
moisture  content  of  the  top  2  cm  with  ground-based  passive  microwave 
data.  The  radiometric  wavelength  bands  tested  were  2.8,  21.4,  and  a  33 
to  50  cm  variable  UHF  band.  They  concluded  that  L-band  sensors  would 
give  the  best  indication  of  volumetric  soil  moisture  in  the  top  2  cm, 
while  UHF  sensors  would  give  the  best  indication  in  the  top  4  cm. 

Satellite  Passive  Microwave  Investigations 
Satellite  Sensors 

The  first  satellite  that  carried  a  passive  microwave  sensor  useful 
for  soil  moisture  experiments  was  launched  in  1972.  Since  then,  four 
additional  satellites  have  been  launched,  with  a  future  planned  launch 
in  mid-1985.  Table  1  gives  a  chronology  of  satellite  passive  microwave 
radiomet ry. 

Nimbus  5,  launched  in  Dec  1972,  was  placed  in  a  circular,  sun-syn¬ 
chronous  orbit  at  an  altitude  of  1112  km.  The  passive  microwave 
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instrument  aboard  was  the  Electrically  Scanning  Microwave  Radiometer 
(ESMR).  The  sensor  scanned  in  a  cross-track  direction  with  a  viewing 
angle  ranging  from  nadir  to  50  degrees.  The  radiometer  had  a  single 
wavelength  of  1.55  cm  and  a  resolution  element  of  25  km  at  nadir.  Its 
primary  objective  was  to  map  surface  features  of  the  earth,  particular¬ 
ly  the  distribution  of  polar  ice  and  precipitating  clouds  over  the 
ocean.  The  ESMR  sensor  provided  about  three  years  of  passive  microwave 
data. 

Skylab  was  the  next  in  the  series,  launched  in  May  1973.  This  was 
the  first  sensor  to  have  a  primary  objective  of  sensing  soil  moisture. 
Skylab  had  a  circular  orbit  at  an  altitude  of  438  km.  The  satellite 
carried  two  passive  microwave  instruments  useful  in  soil  moisture 
experiments.  The  S-193  instrument  mechanically  scanned  in  a  parabola 
at  a  wavelength  of  2.16  cm.  The  smallest  resolution  element  was  16 
km.  The  S-194  instrument  was  a  non-scanning  sensor  with  a  wavelength 
of  21.4  cm.  The  smallest  resolution  element  of  this  L-band  sensor,  115 
km,  was  considerably  larger  than  the  smallest  resolution  element  of 
the  S-193  instrument.  All  soil  moisture  experiments  aboard  Skylab  were 
accomplished  between  May  1973  and  February  1974. 

Nimbus  6  was  launched  in  June  1975  with  an  ESW1  sensor  aboard.  It 
had  a  circular  sun-synchronous  orbit  at  an  altitude  of  1100  km.  This 
sensor  differed  from  the  Nimbus  5  ESMR  in  that  it  had  a  0.81  cm  dual 
polarized  sensor  which  had  a  conical  scan  at  a  constant  angle  of  45 
degrees  with  respect  to  the  earth's  surface.  The  change  in  wavelength 
from  ESMR-5  to  ESMR-6  effectively  tripled  the  instrument's  sensitivity 
to  water  droplets  while  retaining  essentially  the  same  sensitivity  to 
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water  vapor.  Most  of  the  studies  involving  ESMR-6  dealt  with  parame¬ 
ters  over  the  ocean,  such  as  rainfall  rate  or  wind  speed.  Unfortunate¬ 
ly  the  sensor  was  operational  only  for  approximately  two  years. 

Both  Seasat  and  Nimbus  7  were  launched  in  1978.  The  sensor  aboard 
Seasat  obtained  data  between  June  1978  and  October  1978.  It  was  placed 
in  a  circular  near-polar  orbit  at  an  altitude  of  800  km.  Nimbus  7 
replaced  Seasat  in  Oct  1978  and  is  currently  operational  in  a  circular 
sun-synchronous  orbit  at  an  altitude  of  955  km.  Each  satellite  carried 
a  dual  polarized  Scanning  Multichannel  Microwave  Radiometer  (SMMR) 
which  operated  in  five  wavelengths,  0.81,  1.36,  1.66,  2.80,  and  4.54 

cm.  The  sensor  scanned  in  a  conical  pattern  at  a  nearly  constant  angle 
of  about  50  degrees  to  the  earth's  surface.  One  of  the  principal 
parameters  to  be  measured  or  inferred  by  the  SMMR  sensor  was  soil 

moisture. 

The  Defense  Meteorological  Satellite  Program  (DMSP)  satellite  with 
a  planned  launch  in  mid-1985  will  have  a  passive  microwave  radiometer 
aboard.  The  satellite  will  be  placed  in  a  near-polar  sun-synchronous 
orbit  at  an  altitude  of  833  km.  This  DMSP  satellite  will  carry  a 
Special  Sensor  Microwave  Imager  (SSM/I)  described  in  detail  by 
Hollinger  and  Lo  (1983).  This  sensor  will  have  a  conical  scan  at  a 
constant  angle  of  53.1  degrees  to  the  earth's  surface  and  will  operate 

at  wavelengths  of  0.36,  0.81,  1.35,  and  1.55  cm.  All  four  wavelengths 

will  have  a  horizontally  polarized  brightness  temperature  component, 
and  all  but  the  1.35  cm  band  will  have  a  vertically  polarized 
component.  The  resolution  will  range  from  about  15  km  for  the  0.36  cm 
channel  to  about  55  km  for  the  1.55  cm  channel. 
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Resolution 

Several  studies  have  indicated  that  the  lower  resolution  data  from 
satellite-based  sensors  is  sufficient  to  infer  soil  moisture  amounts. 
Hardy  et_  al_.  (1981)  noted  that  operational  soil  moisture  models  used 
in  government  agencies  used  meteorological  data  which  had  a  resolution 
on  the  order  of  100  km.  Newton  et_  ^1_.  (1982)  simulated  the  effect  of 
degrading  the  resolution  from  a  passive  microwave  sensor  for  X,  C,  and 
L-band  wavelengths.  They  found  the  sensitivity  to  soil  moisture 
actually  increased  as  the  simulated  resolution  was  degraded.  The 
optimum  resolution  was  found  to  be  20  km  with  the  sensitivity  of  a  60 
km  resolution  cell  being  equal  to  that  of  a  5  km  cell.  McFarland  and 
Harder  (1982)  simulated  a  degradation  of  resolution  of  Nimbus  5  ESMR 
data  from  25  km  to  50  km.  They  found  the  50  km  correlations  between 
emissivity  and  API  were  as  good  as  or  better  than  the  original  25  km 
resolution  data. 

Statistical  Analysis  Techniques 

The  most  common  method  found  in  the  literature  to  analyze  passive 
microwave  data  statistically  for  soil  moisture  information  was  simple 
linear  regression.  Several  researchers  extended  their  analysis  of 
simply  correlating  the  microwave  data  with  a  measure  of  soil  moisture 
by  using  additional  techniques  such  as  nultiple  linear  regression  or 
stepwise  linear  regression.  Simple  linear  regression  was  used  to 
analyze  the  data  in  many  of  the  Skylab  and  Nimbus  5  experiments,  while 
studies  involving  Nimbus  7  data  tended  to  include  additional 
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techniques.  This  current  trend  has  been  assisted  by  technological 
advances  such  as  improved  parameteri zation  of  the  variables  involved  in 
the  analysis  scheme  resulting  in  reduced  overall  error,  the  advent  of 
higher  resolution  dual  polarization  multi-frequency  microwave  data,  and 
the  increasing  availability  of  advanced  computerized  statistical 
packages. 

Rodgers  and  Si ddal i ngaiah  (1983)  initially  analyzed  three  wave¬ 
lengths  of  dual  polarized  Nimbus  7  SMMR  data  using  simple  linear 
regression  to  delineate  rainfall  areas  over  land.  They  extended  the 
statistical  analysis  with  a  Fisher  linear  discriminant  classification. 
This  was  analagous  to  their  earlier  work  where  they  used  Nimbus  6  ESMR 
data  to  determine  rainfall  over  land  (Rodgers  et_^]_.,  1979). 

Pandey  and  Kakar  (1983)  used  a  multiple  linear  regression  model 
for  retrieval  of  geophysical  parameters  from  SMMR  data.  Spencer  et 
al .  (1983)  applied  a  step-wise  multiple  linear  regression  approach 

which  used  Nimbus  7  SMMR  data  to  infer  rainfall  rates  over  land.  Burke 
and  Ho  (1981)  used  the  Statistical  Parameter  Inversion  Method  (SPIM) 
to  infer  soil  moisture  from  dual  polarized  2.8  and  21  cm  microwave 
data.  The  SPIM  algorithm,  according  to  Hollinger  and  Lo  (1983),  will 
be  used  in  the  analysis  of  future  DMSP  SSM/I  data. 

Investigative  Results 

The  first  available  satellite  sensor  to  be  used  in  investigations 
to  determine  the  relationship  between  satellite  passive  microwave 
brightness  temperatures  and  soil  moisture  was  the  1972  Nimbus  5-ESMR 
sensor.  The  launch  of  Nimbus  6  marked  the  first  time  a  dual  polarized 
sensor  could  be  used  to  determine  this  relationship,  however  the  early 
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demise  of  the  ESMR  sensor  limited  its  investigative  usefulness.  The 
majority  of  investigations  into  this  relationship  currently  found  in 
the  literature  involved  single  polarized  Skylab  data  and  dual  polarized 
Nimbus  5-ESMR  data. 

Stucky  (1975)  used  simple  linear  regression  to  determine  the  cor¬ 
relation  between  Skylab  S-193  brightness  temperatures  and  a  10-day 
API.  This  API  model  used  a  constant  depletion  coefficient  of  0.9.  He 
reported  a  correlation  of  -0.9  which  decreased  nearly  linearly  to  -0.2 
as  the  viewing  angle  was  increased  to  40  degrees.  He  found  also  that 
correlations  for  viewing  angles  less  than  25  degrees  were  higher  for  a 
10-day  API  than  a  6-day  model. 

McFarland  (1976)  used  simple  linear  regression  to  determine  the 
relationship  between  Skylab  S-194  brightness  temperatures  and  an  11-day 
API.  He  reported  a  correlation  of  -0.97  for  June  1973  data  and  -0.73 
for  August  1973  data.  He  suggested  the  decrease  in  the  correlation  was 
due  to  increased  surface  roughness.  Eagleman  and  Lin  (1976)  related 
Skylab  S-194  data  with  actual  soil  moisture  field  samples  collected 
near  the  satellite  overpass  times  with  correlations  similar  to  those  of 
McFarland  (1976). 

Schmugge  et  al_.  (1977)  reported  a  linear  correlation  between 
Nimbus  5-ESMR  data  and  antecedent  rainfall  totals  over  bare  agricultur¬ 
al  soils.  McFarland  and  Blanchard  (1977),  Theis  (1979),  Blanchard  et 
al .  (1981),  and  McFarland  and  Harder  (1982,  1983)  reported  linear  cor¬ 

relations  between  a  normalized  brightness  temperature  and  API.  The 
normalized  brightness  temperature  was  computed  by  dividing  the  bright¬ 
ness  temperature  by  the  daily  maximun  air  temperature.  They  called 
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this  result  an  emissivity  or  a  normalized  emissivity.  The  highest  cor¬ 
relation  coefficients  were  found  with  essentially  bare  agricultural 
soils.  McFarland  and  Harder  (1983)  noted  that  emissivities  over  frozen 
or  snow-covered  soil  were  distinctly  different  than  those  over  bare 
unfrozen  soil.  They  also  suggested  that  soil  moisture  conditions  could 
be  inferred  from  a  stressed  crop  canopy  since  the  vegetative  response 
was  related  to  the  underlying  soil  moisture  conditions.  Harder  (1984) 
found  that  correlations  using  combinations  of  Nimbus  7  SMMR  polarized 
brightness  temperatures  and  API  gave  similar  results  to  correlations 
between  emissivities  and  API.  He  reported  that  seasonal  segments  of 
either  the  emissivity  or  various  transformations  of  polarized  bright¬ 
ness  temperatures  had  higher  correlation  coefficients  with  API  than 
annual  models,  however  annual  models  were  more  reliable  due  to  a  larger 
sample  size. 
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CHAPTER  III 

MATERIALS  AND  METHODS 

Data 


Primary 

The  primary  data  used  in  this  investigation  involved  passive 
microwave  data  from  the  Nimbus  7  Scanning  Multichannel  Microwave  Radio¬ 
meter  (SMMR),  climatic  data  from  the  National  Climatic  Data  Center,  and 
computed  Antecedent  Precipitation  Index  (API)  values  which  were  used  as 
ground  observations. 

NASA's  Goddard  Space  Flight  Center  provided  computer  tapes  with 
Nimbus  7  SMMR  data  from  25  Oct  1978  to  10  Nov  1979.  The  orbit  of  Nim¬ 
bus  7  provided  repetitive  coverage  every  2  or  3  days  over  the  selected 
study  area.  The  local  observation  time  varied  by  about  an  hour  with 
average  overpass  times  of  about  0100L  and  1100L  (0600Z  and  1600Z),  for 
night  and  day  passes  respectively.  The  data  tapes  included  the  nominal 
latitude  and  longitude  of  each  wavelength's  instantaneous  field  of  view 
(IFOV) .  These  coordinates  were  assumed  accurate  and  used  to  align  the 
SMMR  data.  This  investigation  utilized  all  10  available  SMMR  data 
channels,  which  included  both  vertical  and  horizontal  polarized  bright¬ 
ness  temperatures  in  five  wavelengths:  0.81,  1.36,  1.66,  2.80,  and  4.54 
cm. 

Computer  tapes  containing  weather  observations  from  over  500  sta¬ 
tions  in  the  cooperative  observing  network  were  purchased  from  the 
National  Climatic  Data  Center.  Daily  observations  of  maximum  and 


minimum  temperature,  total  precipitation,  snowfall,  and  snowdepth,  were 
extracted  from  these  tapes  for  use  in  the  analysis. 


The  API  was  used  as  ground  observations  in  this  investigation  due 
to  the  lack  of  available  large-scale  soil  moisture  measurements  over 
the  study  area.  Observed  precipitation  values  from  the  climatic  data 
were  used  to  calculate  daily  values  of  the  API. 

Secondary 

Secondary  data  were  used  to  obtain  additional  information  about 
the  study  area.  Images  from  the  geostationary  meteorological  satellite 
(GOES)  were  used  in  delineating  clouds  for  areal  case  study  analysis. 
National  Weather  Service  (NWS)  weather  radar  summary  charts  were  used 
to  locate  precipitating  clouds.  United  States  Geological  Survey  (USGS) 
topographic  and  land-use  maps  along  with  reported  crop  production  sta¬ 
tistics,  were  used  to  select  and  evaluate  grid  cells  for  a  temporal 
case  study  analyses. 


Organization 


Study  Area 

The  study  area  as  shown  in  Fig.  2,  included  most  of  Kansas,  Okla¬ 
homa,  and  the  Panhandle  of  Texas,  along  with  small  portions  of 
Colorado,  New  Mexico,  and  Missouri.  The  base  map  was  generated  locally 
in  a  polar  stereographic  map  projection  true  at  37  degrees  latitude  and 
100  degrees  longitude.  An  analysis  grid,  which  consisted  of  35  rows 
and  40  columns  of  20  x  20  kilometer  cells,  was  overlaid  onto  the  study 
area,  with  each  county  averaging  between  2  and  4  grid  cells.  The  cells 


were  numbered  starting  from  the  top  left-hand  corner  cell.  Areal 
analysis  included  the  complete  study  area,  while  temporal  analysis  was 
performed  on  selected  case  study  grid  cells. 

The  case  study  grid  cells  with  the  approximate  crop  percentage  of 
the  county  selected  for  this  investigation  are  listed  in  Table  2.  The 
crop  categories  listed  in  this  table  are  the  percentage  of  the  total 
county  acreage  that  was  planted  in  the  specified  non- i rri gated  crop. 
These  percentages  were  computed  from  the  county  size  and  crop  acreage 
statistics  reported  in  the  county  agricultural  statistics,  published 
annually  by  the  State  Crop  and  Livestock  Reporting  Services.  The  case 
study  grid  cells  were  chosen  on  the  basis  of  a  predominantly  agricul¬ 
tural  land-use  with  a  high  percentage  of  non-i rrigated  acreage  planted 
in  winter  wheat. 


TABLE  2 

CASE  STUDY  GRID  CELLS  WITH  APPROXIMATE  CROP  PERCENTAGE  OF  COUNTY 


TOTAL 

WHEAT  SORGHUM  CORN  COTTON  HAY  IRRIGATED 


ROW 

COLUMN 

COUNTY 

STATE 

(%) 

(%) 

(%) 

(%) 

(%)  ACREAGE (%) 

03 

27 

CLOUD 

KS 

36 

15 

4 

0 

5 

2 

09 

21 

RUSH 

KS 

36 

8 

1 

0 

4 

2 

09 

29 

MCPHERSON 

KS 

45 

18 

2 

0 

4 

4 

15 

29 

SUMNER 

KS 

66 

5 

0 

0 

3 

0 

18 

27 

GRANT 

OK 

65 

21 

0 

0 

3 

0 

21 

14 

OCHILTREE 

TX 

42 

9 

0 

0 

0 

17 

29 

22 

KIOWA 

OK 

41 

2 

0 

9 

3 

1 

Crop  categories  not  listed  comprised  <  1.5%  of  the  total  county 
acreage. 
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Figure  3  lists  the  1980  dryland  winter  wheat  acreage  by  county  for 
the  study  area.  The  grid  cells  surrounding  the  selected  cells  were 
also  examined  using  this  criteria,  since  the  SMMR  IFOV  of  the  longer 
wavelengths  covered  more  than  one  grid  cell.  A  minimum  of  2  grid  cells 
between  selected  case  study  grid  cells  was  necessary  due  to  the  size  of 
the  microwave  footprint. 

Primary  Database 

The  primary  database  consisted  of  SMMR,  climatic,  and  API  values. 
The  SMMR  and  climatic  data  were  analyzed  objectively  to  the  35  x  40 
grid  using  an  exponential  weighting  function  similar  to  that  used  by 
McFarland  and  Blanchard  (1977).  All  observations  within  a  specified 
radius-of- influence  surrounding  the  grid  cell  were  assigned  weights 
based  on  distance  in  order  to  assign  a  value  at  the  center  of  the 
cell.  The  radius-of-inf luence  for  brightness  temperature,  maximum  and 
minimum  air  temperature,  and  precipitation  were  selected  on  the  basis 
of  density  of  reporting  stations  or  available  SMIW  data.  This  analysis 
algorithm  has  been  reported  in  detail  i  'arder  and  McFarland  (1984). 

The  API  used  in  this  investigation  was  computed  using  an  algorithm 
from  McFarland  and  Harder  (1982).  This  form  was  defined  as: 

API  (i )  =  API (1-1)  *  k  +  P(i )  (12) 

where  k  =  depletion  coefficient 
P  =  effective  precipitation 

The  depletion  coefficient  was  selected  in  a  manner  consistent  with 
McFarland  and  Blanchard  (1977)  and  Harder  (1984).  It  was  defined  as  a 
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Figure  3. 


Percentage  of  1980  county  acreage  planted  In 
dryland  winter  wheat  for  Kansas,  Oklahoma,  and 
Northwest  Texas  (adapted  from  Harder,  1984). 
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sinusoidal  curve  with  a  maximum  value  of  0.92  on  1  February  and  a  mini- 
mun  value  of  0.70  on  1  August.  This  forces  the  depletion  coefficient 

to  be  a  function  of  season,  which  is  related  to  potential  evapotrans- 

piration.  The  effective  precipitation  was  computed  by  raising  the  pre¬ 
cipitation,  in  millimeters,  to  a  power  of  0.891.  The  exponent  was 
initially  suggested  by  Blanchard  et_  a_l_.  (1980)  in  a  study  involving 

runoff  for  watersheds  near  Chickasha,  Oklahoma. 

Snowdepth  was  included  in  the  climatic  database  to  account  for 
accumulated  precipitation  in  the  form  of  snow  which  would  not  immedi¬ 
ately  contribute  to  the  soil  moisture  content.  Once  a  positive  snow- 
depth  had  been  reported,  the  API  was  held  constant,  and  all  subsequent 
precipitation  was  assumed  to  occur  as  snow  until  the  reported  snowdepth 
reached  zero.  Once  the  reported  snowdepth  reached  zero,  all  accumu¬ 
lated  precipitation  was  added  to  the  API  computation  as  effective 

precipitation. 

The  10  channels  of  interpolated  SMW1  brightness  temperatures  were 
placed  in  SMMR  map  files,  one  file  for  each  day's  data.  Each  file  con¬ 
tained  350  records,  one  record  for  each  of  the  35  rows  in  each  of  the 
10  channels  of  brightness  temperature  data.  Each  record  represented 
one  row,  which  contained  40  bytes  of  data.  The  record  length  was  set 
at  40  bytes  since  there  were  40  brightness  temperature  values,  one  per 
grid  cell,  in  each  row.  The  brightness  temperatures  were  expressed  in 
degrees  Kelvin  minus  180,  and  represented  an  integer  between  -128  and 
127,  with  -128  reserved  for  missing  data. 

The  interpolated  climatic  data  values  were  placed  in  climatic  map 
files  similar  to  the  SMMR  map  files,  with  one  file  per  day.  Each  file 


40 


contained  5  channels  of  data:  maximum  air  temperature,  minimum  air 
temperature,  precipitation,  snowfall,  and  snowdepth.  Like  the  SMMR 
files,  each  byte  contained  one  integer  value.  The  climatic  data  had  to 
be  converted  to  a  form  compatible  with  the  computer's  numeric  range  of 
-128  to  127.  Therefore,  the  temperature  data  were  expressed  in  degrees 
Celsius  times  2.5,  precipitation  as  millimeters  minus  126,  and  both 
snowfall  and  snowdepth  as  centimeters  minus  126.  Once  again,  -128  was 
reserved  for  missing  values. 

The  calculated  API  values  between  1  July  1978  and  10  November  1979 
were  placed  in  daily  API  map  files  similar  to  the  SIW  and  climatic  map 
files. 

Image  processor  files  were  created  from  SMMR,  climatic,  API,  and 
GOES  data.  The  GOES  image  was  initially  digitized  to  create  numeric 
intensity  values.  The  numeric  values  from  these  four  data  sources  were 
assigned  an  intensity  level  by  the  image  processor,  and  displayed  as  a 
visual  image.  The  numeric  values  of  the  input  data  sources  could  be 
manipulated  in  a  variety  of  ways  to  display  combinations  or  transforma¬ 
tions  of  the  input  data  to  emphasize  qualitatively  a  certain  feature. 

Secondary  Database 

The  secondary  data  consisted  of  geostationary  meteorological 
satellite  (GOES)  images.  National  Weather  Service  (NWS)  weather  radar 
summary  charts.  United  States  Geological  Survey  (USGS)  topographic  and 
land-use  maps,  and  reported  crop  production  statistics.  The  first  two 
sources  of  secondary  data  gave  an  indication  of  the  amount  of  moisture 
present  in  the  intervening  atmosphere  between  the  object  sensed  and  the 
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radiometer,  which  would  attenuate  the  brightness  temperature  signal. 
The  last  two  data  sources  gave  an  indication  of  the  surface  roughness 
and  vegetation  cover. 


Analysis  Methods 


Temporal  Analysis 

A  computer  program  was  written  to  prepare  individual  grid-cell 
time  series  files  of  the  daily  10  SMMR  brightness  temperature  values, 
daily  values  of  the  5  climatic  data,  and  the  computed  daily  API 
values.  The  program  was  structured  to  compute  and  print  normalized 
brightness  temperature  (emissivity)  values  in  lieu  of  brightness  tem¬ 
perature  values.  Normalized  brightness  temperature  (emissivity)  for 
night  coverage,  was  defined  as  the  ratio  of  the  brightness  temperature 
to  the  minimum  air  temperature.  For  day  passes,  it  was  the  ratio  of 
the  brightness  temperature  to  the  maximun  air  temperature.  This  nor¬ 
malization  of  the  brightness  temperature  has  the  major  function  of 
removing  seasonal  and  shorter  period  temperature  trends  from  the  data, 
but  also  was  assumed  to  approximate  the  emitting- layer  temperature, 
following  work  by  Theis  (1979).  This  program  would  also  plot  any 
series  of  data  for  a  particular  grid  cell  in  order  to  identify  monthly, 
bi-monthly,  or  annual  trends. 

Another  computer  program  was  written  to  analyze  statistically  any 
combination  of  two  of  these  16  available  channels,  for  any  specified 
period  of  time.  The  methodology  used  involved  simple  and  multiple 
linear  regression,  following  previous  work  by  Rodgers  and  Siddalingaiah 
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(1983)  and  Harder  (1984). 

Harder  (1984)  investigated  the  relationship  between  API  and  the 
response  from  54  combinations  of  1.66  cm  wavelength  passive  microwave 
brightness  temperature,  surface  temperature,  and/or  a  sinusoidal 
function  of  the  day  of  the  year.  He  found  that  three  of  these  indices 
accounted  for  most  of  the  variance  in  the  linear  regression  model. 

A  computer  program  was  written  to  use  data  from  any  of  the  5  SMMR 
wavelengths  to  compute  these  three  indices:  horizontally  polarized 
normalized  brightness  temperature  (Th/T ) ,  polarization  difference 
(Tv-Th),  and  normalized  polarization  difference  [Tv-Th/0.5(Tv+Th)]. 
The  program  would  also  compute  statistics  using  these  indices  on  either 
a  monthly,  bi-monthly,  or  annual  basis.  The  bi-monthly  analysis  could 
be  accomplished  by  using  either  the  Jan-Feb  series  Harder  used,  or  a 
Feb-Mar  series,  which  was  thought  to  be  a  better  representation  for  the 
winter  wheat  calendar  during  the  fall  months,  over  this  study  area. 

Areal  Analysis 

The  areal  analysis  utilized  numeric  map  files  spanning  the  entire 
study  area  for  any  one  particular  day.  These  map  files  contained  all 
10  channels  of  SMMR  brightness  temperature  data,  the  5-channel  climatic 
data,  and  the  1-channel  computed  API.  An  example  of  these  files  for 
one  selected  case  study  day,  10  May  1979,  is  given  in  Appendix  A, 
Figures  A-l  through  A-16.  These  numeric  files  were  processed  into 
visual  Images  which  were  used  for  qualitative  comparisons  between  data 
sources  and  selected  channels  within  data  sources.  These  image  proces¬ 
sor  files  were  combined  numerically  in  a  variety  of  ways  to  Illustrate 
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and  help  identify  patterns  not  readily  discernible  on  the  initial 
numeric  map  files. 

The  numeric  map  files,  image  processor  files,  GOES  imagery,  and 
NWS  radar  summary  charts,  were  used  in  the  selection  of  specific  case 
study  days.  Once  these  days  were  chosen,  the  LfSGS  topographic  and 
land-use  maps,  and  the  reported  crop  statistics,  provided  additional 


CHAPTER  IV 


ANALYSIS  AND  DISCUSSION 
Qualitative  Analysis  Methods 


Temporal  Analysis 

The  initial  analysis  consisted  of  a  qualitative  screening  of  a 
grid  cell  16-channel  time  series  file.  The  first  step  in  the  selection 
process  to  determine  which  grid  cell  to  select,  involved  the  use  of  the 
USGS  topographic  and  land-use  maps  for  the  location  of  a  relatively 
flat  agricultural  area  near  the  center  of  the  study  area.  The  reported 
crop  and  county  statistics  from  the  State  Crop  and  Livestock  Reporting 
Service  were  then  examined  to  identify  a  county  with  a  large  percentage 
of  winter  wheat  acreage.  Sumner  County,  Kansas,  located  near  the 
center  of  the  study  area,  was  chosen  since  it  had  the  highest  reported 
county  percentage  of  winter  wheat  acreage  of  the  entire  study  area. 
The  grid  cell  selected,  row  15,  column  29,  was  located  in  the  center  of 
the  northern  half  of  Sumner  County.  Similar  criteria  were  used  to 
select  six  more  agricultural  case  study  grid  cells  spread  across  the 
entire  study  area. 

Areal  Analysis 

The  16-channel  time  series  files  extending  from  1  July  1978  to  10 
Nov  1979  for  these  case  study  grid  cells  were  examined  to  Identity 
Individual  case  study  days  to  be  used  in  an  areal  analysis.  Selection 
criteria  consisted  of  5  or  more  consecutive  days  with  no  reported 
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precipitation,  so  the  observed  microwave  response  would  not  be  unduly 
biased  from  a  recent  rainfall.  The  location  and  amount  of  precipita¬ 
tion  were  taken  from  the  reported  climatic  data,  and  through  analysis 
of  the  NWS  weather  radar  summary  charts.  The  existance  of  clouds  over 
the  study  area  was  determined  through  an  examination  of  the  GOES 
imagery  taken  near  the  SMMR  overpass  times.  Case  study  days  selected 
for  analysis  consisted  of:  days  with  no  precipitation,  days  with  isola¬ 
ted  convective  cells  resulting  in  patchy  areas  of  precipitation,  and 
one  day  with  overcast  conditions  including  a  line  of  thunderstorm  cells 
which  covered  approximately  half  the  study  area.  This  day,  10  May 
1979,  was  selected  as  the  dominant  case  study  day  for  an  areal  analysis 
due  to  the  diverse  atmospheric  conditions. 

Images  of  all  10  SMMR  brightness  temperature  channels  for  10  May 
1979  were  created  through  an  image  processor.  Photographs  were  taken 
of  these  images  to  provide  a  qualitative  comparison  of  the  actual 
brightness  temperature  response  in  these  various  wavelengths  and  polar¬ 
izations.  The  brightness  temperature  ranged  from  approximately  230  K 
to  near  270  K  between  the  10  channels,  with  the  shorter  wavelengths 
recording  the  coldest  temperatures.  The  horizontal  channels  consis¬ 
tently  exhibited  a  larger  temperature  range  for  a  particular  wave¬ 
length,  which  allowed  better  pattern  recognition.  The  varying  resolu¬ 
tions  of  the  SMM*  wavelengths  were  obvious  on  the  images  ranging  from 
high  resolution  on  the  shorter  wavelength  imagery  to  rather  low  resolu¬ 
tion  of  the  longer  wavelength  imagery.  Each  individual  20  x  20  km  grid 
cell  was  evident  on  the  0.81  cm  wavelength  imagery,  while  up  to  12  grid 
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cells  indicated  nearly  identical  temperatures  on  the  4.54  cm  wavelength 
imagery. 

Figure  4a  is  a  reproduction  of  the  GOES  imagery  acquired  near  the 
SMMR  overpass  time.  This  image  shows  the  cloud  cover  present  when  the 
brightness  temperature  data  was  collected,  with  circular  thunderstorm 
cells  indicated  along  the  squall  line.  Figure  4b  is  an  MB  or  thunder¬ 
storm  enhancement  of  the  GOES  image,  which  indicates  the  presence  of 
thunderstorms  over  most  of  the  western  half  of  the  study  area.  The  NWS 
weather  radar  summary  chart  depicted  in  Fig.  5,  indicated  a  squall 
line  with  associated  areas  of  heavy  precipitation  located  northeast- 
southwest,  just  east  of  the  center  of  the  study  area.  The  area  of 
maximum  precipitation  indicated  by  the  climatic  data,  was  located  along 
the  squall  line,  with  a  secondary  maximum  running  from  the  southwestern 
corner  of  Kansas,  northeast  to  the  center  of  the  northern  border  of 
Kansas.  The  recorded  precipitation  amount  was  greatest  near  the  top  of 
the  study  area,  near  the  Kansas  border.  These  two  isolated  areas  of 
heavy  precipitation  made  this  day  a  rather  unique  and  very  informative 
case  study  day. 

Images  created  from  the  actual  numeric  brightness  temperature 
values  greatly  enhanced  the  visual  patterns  indicated  by  both  polariza¬ 
tions  of  the  five  wavelengths  of  the  SMMR  brightness  temperatures. 
Specific  patterns  were  evident  on  each  of  the  different  wavelength  map 
files  of  the  SMMR  brightness  temperatures.  An  explanation  of  the  iden¬ 
tifiable  patterns  seen  on  the  various  wavelength  files  can  be 
provided  by  Fig.  6,  which  gives  the  atmospheric  water  vapor  and  o*y9®n 
attenuation  curves.  The  actual  SMMR  wavelengths  are  given  to  assist  in 
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Figure  5.  NWS  weather  radar  summary  chart  valid  06352, 
10  May  1979. 


FREQUENCY/WAVELENGTH 


Attenuation 
In  dB/km 


Figure  6.  Atmospheric  attenuation  of  microwave  radiation  from  water  vapor 
and  oxygen  (adapted  from  llntz  and  Slmonett,  1976). 
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evaluating  the  influence  of  the  intervening  atmosphere  on  the  response 
from  each  specific  wavelength.  A  comparison  of  both  polarizations  of 
the  0.81  and  1.66  cm  wavelength  imagery  is  provided  by  Fig.  7.  The 
brightness  temperature  pattern  seen  on  both  the  horizontal  and  vertical 
polarization  photographs  of  each  wavelength  illustrates  the  varying 
response  of  the  brightness  temperature  resulting  from  the  effects  of 
the  intervening  atmosphere.  This  pattern  is  more  pronounced  on  the 
horizontal  polarization  photographs,  illustrating  the  opportunity  to 
extract  more  information  from  horizontally  polarized  brightness  temper¬ 
ature  data  in  a  qualitative  soil  moisture  analysis. 

The  patterns  on  both  polarizations  of  the  2.80  and  4.54  cm  wave¬ 
length  imagery  illustrated  by  Fig.  8,  show  the  coldest  temperatures 
near  the  top  center  of  the  study  area  close  to  the  Kansas  border.  This 
was  near  the  area  of  the  maximum  recorded  precipitation  for  the  second¬ 
ary  area,  and  most  probably  the  area  of  overall  maximum  surface  rain¬ 
fall  at  the  time  of  the  SMMR  overpass.  Since  this  pattern  becomes  more 
evident  at  the  longer  wavelengths,  an  explanation  based  on  theory 
involves  the  observation  that  the  longer  wavelengths  are  least  affected 
by  atmospheric  water  vapor,  while  the  effect  of  the  intervening  atmo¬ 
sphere  remains  clearly  evident  of  the  1.66  cm  and  shorter  wavelengths. 
This  hypothesis  suggests  that  the  optimum  wavelength  to  use  to  sense 
surface  precipitation  or  surface  soil  moisture  without  significant  loss 
of  resolution,  would  be  the  2.80  cm  wavelength.  While  the  patterns  on 
the  2.80  and  4.54  cm  imagery  were  very  similar,  the  degraded  resolution 
of  the  4.54  cm  imagery  was  quite  pronounced.  This  significantly 


Figure  7-  SMMR  brightness  temperature  images  for  0.81  and  1.66  cm  horizontal  and  vertical 


Figure  8.  SMMR  brightness  temperature  images  for  2.80  and  4.54  cm  horizontal  and  vertical 
polarizations  for  10  May  1979  for  the  study  area. 
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reduced  resolution  emphasizes  the  usefulness  of  the  shorter  wavelength 
imagery  when  examining  large  scale  patterns. 

The  pattern  on  the  0.81  cm  wavelength  imagery  shown  on  Figs.  9 
and  10,  identified  the  area  of  maximum  radar  backscatter  along  the 
squall  line,  which  indicated  the  maximum  area  of  precipitating  clouds. 
A  secondary  nearly  circular  area  in  southwestern  Kansas  was  also 
evident  on  the  imagery.  This  area  was  covered  by  precipitating  clouds, 
but  received  approximately  one- fourth  the  total  rainfall  amount  of  the 
area  near  the  squall  line.  A  likely  explanation  of  the  observed  phen¬ 
omena  is  that  the  0.81  cm  wavelength  signal  was  strongly  attenuated  in 
both  cases  by  the  precipitating  clouds.  The  area  of  maximum  radar 
backscatter  along  the  squall  line  was  clearly  discernable  only  on  the 
0.81  cm  wavelength  imagery.  The  pattern  on  the  1.36  cm  wavelength 
imagery  gave  little  indication  of  Tower  temperature  over  the  area  of 
maximum  precipitating  clouds,  but  indicated  the  area  of  secondary  pre¬ 
cipitation  in  southwestern  Kansas.  This  wavelength  is  in  an  atmo¬ 
spheric  water  vapor  attenuation  region  shown  by  Fig.  6,  which  explains 
the  relatively  uniform  pattern  indicated  by  this  imagery.  This  wave¬ 
length  is  of  great  usefulness  in  detecting  the  presence  of  atmospheric 
water  vapor,  but  is  of  limited  usefulness  for  an  investigation  into  the 
relationship  between  passive  microwave  brightness  temperatures  and  API, 
or  soil  moisture.  Since  this  wavelength  had  the  smallest  temperature 
range  of  all  5  wavelengths,  which  caused  the  relatively  uniform  pattern 
across  the  image,  it  was  excluded  from  the  figures  provided  in  this 
study.  The  secondary  area  of  maximum  precipitation  detectable  on  the 
1.36  cm  wavelength  imagery  was  defined  more  vividly  on  the  1.66  cm 
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wavelength  imagery  without  a  noticeable  loss  of  resolution.  A 
suggested  explanation  for  this  observation  is  that  the  brightness 
temperatures  from  the  0.81  cm  imagery  were  strongly  attenuated  by  the 
precipitating  clouds,  while  the  1.36  cm  and  longer  wavelengths  provided 
a  better  indication  of  the  actual  surface  moisture  conditions.  This 
explanation  was  weighted  by  the  fact  that  the  reported  rainfall 
indicated  by  the  climatic  data  was  a  24  hour  total,  while  the  passive 
microwave  imagery  provided  an  instantaneous  response  at  the  time  of  the 
SMMR  overpass. 

Quantitative  Temporal  Analysis 

The  quantitative  temporal  analysis  involved  correlating  various 
transformations  of  the  SMMR  brightness  temperatures  with  the  API.  The 
correlation  coefficients  were  compared  using  the  Fisher's  V 
transformation  described  by  Brooks  and  Carruthers  (1951).  The  equation 
is 

V  =  0.5  In  (1+r/l-r)  (13) 

where  r  is  any  correlation  coefficient.  This  V  statistic  has  a 
nearly  normal  distribution  and  a  standard  error  independent  of  the  V 
value.  The  standard  error  is  given  by 

S.E.Z'  *  1/  fTT  (14) 


where  n  is  the  number  of  observations  involved  In  computing  the 
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correlation  coefficient.  The  calculated  value  of  V  can  be  used  to 
test  the  statistical  significance  of  the  correlation  coefficients.  To 
accomplish  this,  the  S.E.Z'  is  multiplied  by  the  value  of  the  standard 
deviation  for  the  selected  significance  level  (i.e.  1.96  for  a  0.05 
significance  level).  This  result  is  added  to  or  subtracted  from  the  V 
value  to  obtain  the  upper  and  lower  Z'  bounds.  These  V  bounds  can 
then  be  converted  to  correlation  coefficient  bounds  at  a  particular 
significance  level,  by  using  the  inverse  of  equation  (13). 

Verification  of  the  API  Effective  Precipitation  Exponent 

The  effective  precipitation  exponent  used  to  convert  precipitation 
to  effective  precipitation  in  the  calculation  of  the  API  was  0.891. 
This  partitions  rainfall  into  the  soil  moisture  component  (effective 
rainfall)  and  runoff.  The  greater  the  rainfall,  the  greater  the  runoff 
as  a  function  of  the  value  of  the  exponent.  The  significance  of  this 
exponent  for  use  in  this  investigation  was  evaluated  by  varying  the 
exponent  between  0.85  and  0.95  on  a  monthly  basis  for  the  1.66  cm  wave¬ 
length  data  for  row  15  column  29.  These  data  are  presented  in  Appendix 
B,  Table  B-l.  The  monthly  variation  in  the  correlation  coefficient 
between  the  1.66  cm  horizontally  polarized  normalized  brightness  tem¬ 
perature  and  API,  shown  in  Figs.  11a  and  lib,  was  found  to  be  depen¬ 
dent  on  the  amount  and  type  of  the  observed  monthly  precipitation. 
During  the  months  when  the  precipitation  was  very  high,  between  March 
and  May,  the  correlation  coefficient  remained  equally  high  between  all 
four  variations  of  the  effective  precipitation  exponent.  When  the 
rainfall  became  convective  and  decreased  in  total  amount,  between  May 
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Row  15  Col  29  •  *  0.850 


MONTH 


Figure  11a.  Comparison  of  the  1.66  cm  wavelength  correlation  coefficients  between  the 
horizontally  polarized  normalized  brightness  temperature  and  four  API 
effective  precipitation  exponents  for  March  through  July  1979. 
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Figure  lib.  Comparison  of  the  1.66  cm  wavelength  correlation  coefficients  between  the 
horizontally  polarized  normalized  brightness  temperature  and  four  API 
effective  precipitation  exponents  for  August  through  November  1979  and 
all  months  combined. 
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and  August,  the  difference  in  the  correlation  coefficient  between  the 
four  examples  greatly  increased.  Months  with  relatively  high  convec¬ 
tive  rainfall  benefited  from  the  larger  effective  precipitation  expo¬ 
nent,  while  months  with  low  convective  precipitation  benefited  from  the 
lower  exponent.  With  convective  precipitation,  large  amounts  of  rain¬ 
fall  often  occur  over  a  small  area  during  a  short  period  of  time. 
Often,  much  of  this  rainfall  contributes  to  high  rates  of  surface  run¬ 
off,  therefore,  the  larger  effective  precipitation  exponent  would  bet¬ 
ter  account  for  the  occurrance  of  this  type  of  precipitation.  Little 
difference  in  the  correlation  coefficient  was  observed  between  Septem¬ 
ber  and  November  when  precipitation  was  low.  During  light  rainfall 
periods,  little  difference  in  the  way  the  precipitation  was  reduced  to 
effective  precipitation  should  be  evident  since  the  initial  values  are 
minimal.  The  conclusion  is  that  an  effective  precipitation  exponent  of 
0.891  is  a  good  exponent  to  use  if  one  number  had  to  be  chosen  for  the 
entire  year,  although  a  value  of  0.90  appears  to  be  equally  effective. 
A  better  approach  would  be  to  vary  the  exponent  monthly  based  on  a 
climatological  and  topographic  grid  cell  analysis. 

The  correlation  coefficients  for  the  four  effective  precipitation 
exponents  were  compared  using  the  Fisher's  Z'  transformation.  July  and 
August  had  the  greatest  range  in  correlation  coefficients,  however  no 
significant  difference  was  found  between  the  correlation  coefficients 
for  any  of  the  four  effective  precipitation  exponents  for  either  month, 
at  the  0.05  significance  level.  The  differences  in  the  correlation 
coefficients  between  the  months  were  greater  than  within  months,  with 
some  months  such  as  April  and  May  being  significantly  different. 


Correlation  Analysis  on  Transforms 


The  quantitative  temporal  analysis  involved  correlating  three 
transforms  of  the  SMMR  brightness  temperatures  with  the  API.  These 
transforms  were  the  horizontally  polarized  normalized  brightness  tem¬ 
perature  (Th/T),  the  polarized  difference  transform  (Tv-Th),  and  the 
normalized  polarization  difference  transform  [Tv-Th/0.5(Tv+Th)].  The 
data  used  to  calculate  these  correlations  are  given  in  Appendix  B, 
Tables  B-2  to  B-8. 

Figure  12  provides  a  comparison  between  the  average  correlation 
coefficient  between  the  horizontally  polarized  normalized  brightness 
temperature  (Th/T)  and  the  API  for  the  months  of  April  and  September 
1979.  The  correlation  coefficient  is  the  average  of  all  case  study 
grid  cells  except  Ochiltree  County,  Texas,  (row  21,  column  14)  due  to 
the  presence  of  a  large  amount  of  irrigated  acreage.  This  grid  cell 
was  excluded  since  irrigated  soil  would  contaminate  the  SMMR  brightness 
temperature  correlation  with  API  since  the  brightness  temperatures 
would  indicate  a  substantial  recent  rainfall,  which  would  not  be 
evident  on  the  calculated  API.  April  and  September  were  selected  since 
they  represent  months  of  predominantly  bare  soil,  which  is  shown  by  the 
Oklahoma  Crop  Calendar  given  in  Table  3.  The  best  correlation  between 
the  SMMR  brightness  temperature  transforms  and  API  should  occur  under 
these  conditions  since  the  precentage  of  bare  soil  and  intensity  of  the 
precipitation  represent  the  primary  influences  of  the  correlation 
between  the  SMMR  brightness  temperature  and  API.  The  amount  and  inten¬ 
sity  of  precipitation  affect  the  correlation  coefficient  since  the 
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Figure  12.  Comparison  of  the  correlation  coefficients  between  the 
horizontally  polarized  normalized  brightness  temperature 
and  API.  The  correlation  coefficients  are  the  average  of 
six  grid  cells  for  all  SMMR  wavelengths  for  April  (n=57) 
and  September  (n=65). 
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Figure  13.  Comparison  of  the  correlation  coefficients  between  the 

polarization  difference  transform  and  API.  The  correlation 
coefficients  are  the  average  of  six  grid  cells  for  all  SMMR 
wavelengths  for  April  (n-57)  and  September  (n*65). 


OKLAHOMA  CROP  CALENDAR 


Note:  Indicated  periods  relate  to  when  5  to  95%  of  the  crop  is  planted  or  harvested. 
(Adapted  from  Oklahoma  Crop-Weather  Bulletin) 
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sinusoidal  moisture  depletion  is  best  suited  under  conditions  of  high 
amounts  of  uniform  non-convective  rainfall.  Under  dry  conditions,  the 
evapotranspi ration  curve  becomes  more  linear  causing  a  decrease  in  the 
correlation  between  the  API  and  soil  moisture.  The  correlation  coef¬ 
ficients  were  very  high  for  these  months,  ranging  between  0.84  and  0.90 
for  April,  and  between  0.64  and  0.72  for  September.  The  1.66  cm  wave¬ 
length  correlation  coefficient  was  the  highest  among  all  five  wave¬ 
lengths  for  both  months,  but  not  significantly  different  from  the  rest 
of  the  correlation  coefficients,  at  the  0.05  level.  Since  the  effec¬ 
tive  precipitation  exponent  was  developed  using  1.55  cm  wavelength  pas¬ 
sive  microwave  data,  the  1.66  cm  wavelength  would  be  expected  to  have 
the  highest  correlation  with  API  from  the  five  SMMR  wavelengths,  while 
a  different  effective  precipitation  exponent  may  work  better  for 
another  wavelength.  All  wavelengths  of  the  horizontally  polarized 
normalized  brightness  temperature  had  a  correlation  coefficient  near 
0.90  over  bare  ground  for  all  6  case  study  grid  cells. 

Figure  13  gives  the  same  comparisons  as  Fig.  12,  but  with  the 
polarization  difference  (Tv-Th)  transform.  The  correlation  coeffic¬ 
ients  for  April  range  between  0.77  and  0.91,  and  between  0.61  and  0.87 
for  September,  The  correlation  coefficients  from  the  1.36  and  1.66  cm 
wavelength  data  decreased  considerably,  particularly  in  September, 
between  this  transform  and  the  horizontally  polarized  normalized 
brightness  temperature.  The  correlation  coefficients  were  compared 
with  each  other  using  the  Fisher's  V  transformation.  For  the  month  of 
April,  the  2.80  cm  correlation  coefficient  was  the  highest,  but  it  was 
not  significantly  different  from  either  the  0.81  or  4.54  cm  wavelength 


correlation  coefficient,  at  the  0.05  level.  The  1.36  cm  wavelength  had 
the  lowest  correlation  coefficient,  but  not  significantly  different 
from  the  1.66  cm  correlation  coefficient.  For  September,  the  2.80  cm 
wavelength  correlation  coefficient  was  significantly  higher  than  any 
other  coefficient,  while  the  remaining  four  coefficients  did  not  sig¬ 
nificantly  differ  from  each  other.  A  comparison  between  wavelengths  of 
the  polarization  difference  transform  showed  that  the  correlation  coef¬ 
ficients  from  these  wavelengths  were  nearly  identical  to  those  of  the 
normalized  polarization  difference  transform. 

The  correlation  coefficients  for  the  normalized  polarization  dif¬ 
ference  [Tv-Th/0.5(Tv+Th)]  transform,  shown  by  Fig.  14,  ranged  between 
0.73  and  0.91  for  April,  and  between  0.60  and  0.87  for  September.  The 
correlation  coefficients  were  compared  with  each  other  using  the 
Fisher's  V  transformation.  For  the  month  of  April,  the  2.80  cm  cor¬ 
relation  coefficient  was  the  highest,  but  it  was  not  significantly  dif¬ 
ferent  from  either  the  0.81  or  4.54  cm  wavelength  correlation  coef¬ 
ficient,  at  the  0.05  level.  The  1.36  cm  wavelength  had  the  lowest  cor¬ 
relation  coefficient,  but  not  significantly  different  from  either  the 
0.81  or  the  1.66  cm  correlation  coefficient.  For  September,  the  2.80 
cm  wavelength  correlation  coefficient  was  significantly  higher  than  any 
other  coefficient,  while  the  remaining  four  coefficients  did  not 
significantly  differ  from  each  other. 

Figure  15  is  the  correlation  coefficient  from  the  multiple  linear 
regression  on  all  three  transforms.  These  correlation  coefficients 
ranged  between  0.88  and  0.93  for  April,  and  between  0.78  and  0.92  for 
September.  While  the  2,80  cm  wavelength  again  had  the  highest 
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Figure  15.  Comparison  of  the  correlation  coefficients  between  multiple  linear 
regression  on  all  three  transforms  and  API.  The  correlation 
coefficients  are  the  average  of  six  grid  cells  for  all  SMMR 
wavelengths  for  April  (n«57)  and  September  (n-65). 
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correlation  coefficient  for  the  month  of  April,  it  was  not  significant¬ 
ly  higher  than  any  of  the  other  correlation  coefficients,  at  the  0.05 
level.  For  September,  the  2.80  cm  wavelength  correlation  coefficient 
was  significantly  higher  than  all  other  wavelengths,  except  the  4.54 
cm.  No  significant  difference  was  found  between  the  0.81,  1.36,  and 
1.66  cm  wavelength  correlation  coefficients.  While  the  2.80  cm  wave¬ 
length  correlation  coefficient  was  the  highest  for  both  months,  the 
magnitude  of  this  coefficient  for  both  the  polarization  difference  and 
the  normalized  polarization  difference  transforms  was  approximately 
equal.  This  indicates  that  either  one  of  these  transforms  of  the  SMMR 
brightness  temperatures  can  be  used  independently  to  estimate  the  API 
over  bare  soil  without  the  need  of  an  estimate  of  the  surface  tempera¬ 
ture. 

Comparison  of  Transforms  and  Multiple  Linear  Regression 

Figures  16a  and  16b  compare  the  correlation  coefficient  for  each 
of  the  three  transforms  and  the  multiple  linear  regression  on  all  three 
transforms  for  all  five  wavelengths  for  both  April  and  September  re¬ 
spectively.  The  relative  value  of  each  transform  indicated  by  the  mag¬ 
nitude  of  its  correlation  coefficient  is  immediately  apparent,  illus¬ 
trating  the  apparent  superiority  of  the  2.80  cm  wavelength.  An  explan¬ 
ation  for  this  observation  deals  with  the  atmospheric  attenuation  of 
each  wavelength  and  the  degrading  resolution  as  the  wavelength 
increases.  The  0.81  cm  wavelength  is  attenuated  by  both  atmospheric 
water  vapor  and  oxygen,  thereby  correlating  better  with  precipitating 
clouds  rather  than  API  or  soil  moisture.  The  1.36  cm  wavelength  is  in 
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ure  16a.  Comparison  of  the  April  correlation  coefficients  between  the  three 

transforms  and  multiple  linear  regression  (ml r )  on  all  three  transforms, 
and  API.  The  correlation  coefficients  are  the  average  of  six  grid  cells 
{n=S7 )  for  all  SMMR  wavelengths. 
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Figure  16b.  Comparison  of  the  September  correlation  coefficients  between  the  three 
transforms  and  multiple  linear  regression  (mlr)  on  all  three  transforms, 
and  API.  The  correlation  coefficients  are  the  average  of  six  grid  cells 
(n*65)  for  all  SMMR  wavelengths. 
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the  center  of  a  water  attenuation  region  and  therefore  would  be  expect¬ 
ed  to  indicate  the  observed  lower  correlation  with  API.  The  1.66  cm 
wavelength  still  has  a  rather  high  amount  of  atmospheric  water  vapor 
attenuation,  but  slightly  lower  than  the  1.36  cm  band.  Therefore,  the 
correlation  coefficient  between  the  1.66  cm  wavelength  and  API  should 
be,  and  was  found  to  be,  slightly  higher  than  the  1.36  cm  correlation 
coefficient.  Both  the  2.80  and  4.54  cm  wavelengths  are  at  the  lower 
end  of  both  the  atmospheric  water  vapor  and  oxygen  attenuation  curves 
(see  Fig.  6),  and  should  correlate  well  with  API.  The  extremely  high 
correlation  coefficient  from  the  2.80  cm  wavelength  was  consistently 
observed  with  all  case  study  grid  cells  throughout  the  entire  year. 
The  correlation  coefficient  from  the  4.54  cm  wavelength  was  lower  than 
that  derived  when  using  the  2.80  cm  band,  which  is  attributed  to  the 
reduced  resolution  of  the  longer  wavelength,  which  will  dilute  the 
individual  grid  cell  response. 

Comparison  of  Case  Study  Grid  Cells 

The  monthly  correlation  coefficient  between  the  horizontally 
polarized  normalized  brightness  temperature  and  API  for  the  1.66  and 
2.80  cm  wavelengths  for  the  predominantly  wheat  grid  cell,  row  15 
column  29,  is  shown  in  Fig.  17a.  Figures  18a-e  give  the  annual  API 
trace,  which  represents  soil  moisture,  for  this  grid  cell.  The  monthly 
variations  in  the  correlation  coefficient  are  primarilly  the  result  of 
changes  in  amount  ^ f  both  vegetation  and  rainfall.  The  actual  correla¬ 
tion  coefficient  remained  between  0.75  and  0.90  for  the  months  of  March 
and  April  as  a  function  of  the  large  rainfall  and  developing  wheat 
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Figure  17a.  Comparison  of  the  correlation  coefficients  between  the 
horizontally  polarized  normalized  brightness  temperature 
and  API  for  row  15  col  29.  The  correlation  coefficients 
are  for  both  the  1.66  and  2.80  cm  SMMR  wavelengths. 
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.  Comparison  of  the  correlation  coefficients  between  the 
horizontally  polarized  normalized  brightness  temperature 
and  API  for  row  09  col  29.  The  correlation  coefficients 
are  for  both  the  1.66  and  2.80  cm  SMMR  wavelengths. 
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Figure  18c.  Time  series  plots  of  API  for  row  15  col  29 
between  01  Apr  1979  and  20  Jun  1979. 


Figure  18d.  Time  series  plots  of  API  for  row  15  col  29 
between  20  Jun  1979  and  08  Sep  1979. 


Figure  I8e.  Time  series  plots  of  API  for  row  15  col  29 
between  08  Sep  1979  and  10  Nov  1979. 


canopy.  By  May,  the  wheat  had  reached  maturity,  and  the  amount  of  the 


convective  precipitation  was  less  than  either  March  or  April.  The 
large  drop  in  the  correlation  coefficient,  from  >  0.85  to  near  0.40  was 
due  to  the  increasing  reflection  from  the  developing  wheat  canopy.  The 
wheat  harvest  occurred  in  June  and  the  amount  of  rainfall  increased. 
The  increase  in  bare  soil  observed  by  the  SMMR  sensor  contributed  to  a 
large  rise  in  the  correlation  coefficient,  which  increased  to  >  0.90. 
The  low  precipitation  combined  with  the  effects  of  the  miscellaneous 
vegetation  present  in  the  grid  cell,  primarilly  alfalfa  hay  (see  Table 
3),  caused  the  correlation  coefficient  to  remain  near  0.60  for  July  and 
August.  By  October,  as  a  consequence  of  the  continued  dry  weather,  the 
vegetative  cover  declined  and  the  correlation  coefficient  increased  to 
>  0.80.  By  November,  the  amount  of  precipitation  increased  resulting 
in  the  development  of  the  winter  wheat  canopy  and  a  slight  decrease  of 
the  correlation  coefficient  to  approximately  0.75.  The  months  of 
December,  January,  and  February  are  excluded  since  the  observed  bright¬ 
ness  temperatures  are  not  related  to  soil  moisture  when  the  soil  sur¬ 
face  is  either  frozen  or  covered  with  snow. 

The  monthly  correlation  coefficient  between  the  horizontally 
polarized  normalized  brightness  temperature  transform  and  API  for  the 
1.66  and  2.80  cm  wavelengths  for  the  mixed  wheat  and  sorghum  grid  cell, 
row  09  column  29,  is  shown  in  Fig.  17b.  The  associated  annual  API 
trace  is  given  in  Figs.  19a-e.  The  correlation  coefficients  for  March 
and  April  are  extremely  high  (near  0.90)  due  to  the  predominance  of 
bare  soil,  short  vegetation,  and  abundant  rainfall.  By  May,  the  wheat 
had  reached  maturity  and  the  precipitation  decreased,  and  the 


Figure  19a.  Time  series  plots  of  API  for  row  09  col  29 
between  23  Oct  1978  and  12  Jan  1979. 


Figurr  19b.  Time  series  plots  of  API  for  row  09  col  29 
between  12  Jan  1979  and  01  Apr  1979. 


Figure  19c.  Time  series  plots  of  API  for  row  09  col  29 
between  01  Apr  1979  and  20  Jun  1979. 


Figure  19e.  Time  series  plots  of  API  for  row  09  col  29 
between  08  Sep  1979  and  10  Nov  1979. 
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correlation  coefficient  plummetted  to  <  0.10.  In  June,  the  winter 
wheat  was  harvested,  with  the  increase  in  amount  of  bare  soil  seen  by 
the  SMMR  sensor  resulting  in  a  increase  in  the  correlation  coefficient 
to  >  0.90  for  the  1.66  cm  wavelength.  August  had  a  very  low  recorded 
precipitation  amount  and  predominantly  bare  wheat  fields,  however, 
irrigation  of  the  new  seedbed  and  the  maturing  sorghum  crop  caused  the 
correlation  coefficient  to  drop  to  approximately  0.25.  By  September, 
the  correlation  coefficient  had  risen  to  approximately  0.90,  and  stayed 
between  0.80  and  0.90  through  November.  This  increase  was  the  result 
of  the  increased  bare  soil  after  the  summer  crops  were  harvested  and 
the  relatively  poor  condition  of  the  wheat  crop  due  to  the  dry  soil. 
The  correlation  between  the  SMMR  brightness  temperatures  and  the  API  is 
significantly  better  over  bare  soil.  With  knowledge  of  the  surface 
conditions,  the  changes  observed  in  the  correlation  coefficients  can  be 
interpreted  throughout  the  year.  The  correlation  remains  better  over  a 
uniform  agricultural  surface  since  there  are  fewer  interfering  factors, 
although  the  correlation  coefficient  for  a  mixed  agricultural  grid  cell 
still  reaches  the  extremely  high  levels  seen  in  the  predominantly  wheat 
grid  cell.  The  thickness  and  height  of  vegetation  will  strongly  affect 
the  correlation  coefficient  however,  the  type  of  vegetation  in  a  mixed 
agricultural  grid  cell  seems  to  be  irrelevant. 

Supporting  Data 

The  brightness  temperature  correlation  coefficient  and  the  normal¬ 
ized  brightness  temperature  correlation  coefficient  matrices  for  all 
five  wavelengths  and  both  polarizations  for  row  15  column  29,  are  given 
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in  Appendix  C,  Tables  C-l  and  C-2.  The  correlation  coefficient  between 
both  the  horizontally  and  vertically  polarized  brightness  temperatures 
and  API,  for  all  7  case  study  grid  cells,  is  given  in  Appendix  C, 
Tables  C-3  and  C-4.  The  correlation  coefficient  between  the  1.66  cm 
wavelength  transforms  and  API  for  both  a  Jan-Feb  and  Feb-Mar  series  are 
given  in  Appendix  C,  Tables  C-5  and  C-6. 
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CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 
Summary 

The  major  purpose  of  this  research  was  to  find  ways  to  estimate 
the  API  using  nultispectral  dual  polarized  passive  microwave  data. 
Since  the  API  is  by  definition  strongly  related  to  surface  soil 
moisture,  this  investigation  encompasses  estimation  of  surface  soil 
moisture.  This  investigation  involved  selecting  and  qualitatively 
analyzing  case  study  days  spanning  the  entire  study  area.  Photographs 
of  the  recorded  brightness  temperatures  for  a  selected  case  study  day 
were  taken  to  assist  in  a  qualitative  areal  pattern  recognition 
analysis  of  all  10  recorded  channels  of  microwave  data.  A  quantitative 
temporal  analysis  was  accomplished  using  correlation  coefficients 
between  various  brightness  temperature  transforms  and  API  for  the  case 
study  grid  cells  on  a  monthly,  bi-monthly,  and  annual  basis. 

Conclusions 

The  brightness  temperature  range  for  all  five  wavelengths  was 
larger  in  the  horizontal  polarized  channel  than  the  vertical  channel. 
This  allowed  a  more  in-depth  pattern  recognition  analysis.  The  coldest 
brightness  temperature  patterns  corresponded  to  areas  of  maximum  pre¬ 
cipitation,  rather  than  areas  of  maximum  recorded  rainfall. 

The  0.81  cm  wavelength  was  strongly  attenuated  by  precipitating 
clouds  and  correlated  well  with  atmospheric  water  vapor  content.  The 
resulting  pattern  on  the  0.81  cm  wavelength  imagery  indicated  the 
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presence  of  clouds  and  surface  moisture  for  areas  where  clouds  did  not 
attenuate  the  passive  microwave  signal.  The  relatively  small  bright¬ 
ness  temperature  range  from  the  1.36  cm  wavelength  resulted  in  a  very 
uniform  visual  pattern  on  the  brightness  temperature  imagery.  Since 
this  wavelength  is  in  an  atmospheric  water  vapor  attenuation  region, 
1.36  cm  data  would  be  of  value  in  detecting  the  presence  of  atmospheric 
water  vapor,  but  is  of  limited  usefulness  for  a  soil  moisture  inves¬ 
tigation.  The  1.66  cm  wavelength  data  provided  increased  pattern 
detail  and  responded  better  to  the  surface  moisture  conditions  than  the 
shorter  wavelength  data.  The  2.80  cm  wavelength  data  provided  the  most 
distinct  high  resolution  brightness  temperature  pattern  of  all  five 
available  SMMR  wavelengths.  Since  the  horizontal  polarized  channel 
provided  a  larger  brightness  temperature  range,  the  2.80  horizontally 
polarized  channel  would  be  the  best  single  SMMR  channel  to  use  for  a 
surface  soil  moisture  investigation.  The  4.54  cm  wavelength  data 
responded  very  well  to  the  surface  moisture  conditions,  however  the 
lower  resolution  reduced  the  value  of  this  data  for  a  small  area  soil 
moisture  investigation. 

The  API  effective  precipitation  exponent  of  0.891  used  to  convert 
precipitation  to  effective  precipitation  (soil  moisture  component) 
worked  well  for  this  study  area,  although  a  value  of  0.90  had  essen¬ 
tially  the  same  results.  The  topography  and  amount  and  intensity  of 
the  precipitation  seemed  to  be  the  predominant  factors  to  consider  when 
optimizing  an  API  effective  precipitation  exponent  for  a  particular 
area. 

The  correlation  coefficients  between  the  horizontally  polarized 
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normalized  brightness  temperature  (Th/T)  and  API  for  the  months  of 
April  and  September  were  very  high  and  indicated  little  difference 
among  the  five  SMMR  wavelengths.  The  correlation  coefficients  for  this 
transform  over  bare  ground  were  near  0.90  for  all  six  case  study  grid 
cells.  The  correlation  coefficient  between  the  polarization  difference 
(Tv-Th)  transform  and  API  had  a  larger  range  of  values  than  the  pre¬ 
vious  transform,  and  a  slightly  higher  overall  value  among  the  five 
wavelengths  for  both  April  and  September.  The  2.80  cm  wavelength  had 
the  highest  correlation  coefficient  for  all  five  wavelengths.  The  cor¬ 
relation  coefficients  between  the  normalized  polarization  difference 
[Tv-Th/0.5(Tv+Th)]  and  API  were  nearly  identical  to  those  from  the 
polarization  difference  transform,  with  the  2.80  cm  wavelength  again 
having  the  highest  correlation  coefficient  for  all  five  wavelengths  for 
the  months  of  April  and  September.  The  correlation  coefficients  from 
the  multiple  linear  regression  performed  on  all  three  transforms  were 
slightly  higher  than  those  from  any  of  the  individual  transforms.  The 
2.80  cm  wavelength  again  had  the  highest  correlation  coefficient,  how¬ 
ever  this  value  was  approximately  the  same  as  the  value  from  either  the 
polarization  difference  transform  or  the  normalized  polarization  dif¬ 
ference  transform.  This  indicates  that  either  one  of  these  transforms 
could  be  used  independently  to  estimate  the  API  over  bare  soil  without 
an  estimation  of  the  surface  temperature. 

Monthly  changes  in  the  correlation  coefficients  between  the  SMMR 
brightness  temperatures  and  API  can  be  interpreted  throughout  the  year 
using  supplementary  surface  observations.  The  correlation  coefficients 
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are  consistently  higher  over  bare  soil  due  to  fewer  interfering  fac¬ 
tors,  such  as  the  thickness  and  height  of  the  vegetative  cover. 

Recommendations  for  Further  Research 

The  API  model  can  be  improved  upon  by  optimizing  the  effective 
precipitation  exponent  used  in  reducing  precipitation  to  effective  pre¬ 
cipitation.  By  accounting  for  terrain,  type  and  amount  of  vegetation, 
soil  type  and  structure,  and  type  and  rate  of  precipitation,  an 
improved  effective  precipitation  exponent  could  be  developed  on  a  sin¬ 
gle  grid  cell  basis.  An  improvement  over  the  methodology  used  in  this 
investigation  would  be  to  vary  the  exponent  monthly  by  either  grid  cell 
or  latitude.  A  two-layer  model  of  the  effective  precipitation  curve 
could  be  implemented  also  since  the  soil  dry-down  occurs  in  various 
stages  depending  on  the  available  soil  moisture. 

Atmospheric  effects  were  not  specifically  accounted  for  in  this 
investigation.  The  0.81  cm  wavelength  was  shown  to  indicate  the 
presence  of  precipitating  clouds  while  the  1.36  cm  band  indicated 
atmospheric  water  vapor  content.  This  information  could  be  used  to 
modify  the  response  of  the  longer  wavelengths  when  inferring  surface 
soil  moisture. 

A  study  by  Smith  and  Newton  (1983)  used  simulated  passive  micro- 
wave  data  to  drive  a  soil  moisture  model  called  CONSERVB.  With  the 
addition  of  measured  soil  moisture  data,  the  usefulness  of  this  model 
could  be  tested  with  multispectral  dual  polarized  SMMR  data,  to  deter¬ 
mine  whether  the  CONSERVB  or  API  model  correlates  better  with  soil 
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The  SSM/I  sensor  aboard  the  next  DMSP  satellite  Mill  contain 
similar  wavelengths  and  polarizations  as  the  SMMR  sensor.  The  analysis 
methodology  used  in  this  Investigation  will  be  applicable  to  all  future 
investigations  Involving  this  type  of  data.  The  conclusions  drawn  from 
this  investigation  should  be  used  as  a  launching  platform  for  analyzing 
future  passive  microwave  data  utilizing  advanced  statistical 
techniques. 
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APPENDIX  A 

AREAL  MAPS  OF  SMMR,  CLIMATIC,  AND  API  DATA 
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Figure  A-l.  Areal  map  of  the  10  May  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  0.81  cm  horizontally  polarized  wavelength. 
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Figure  A-2.  Areal  map  of  the  10  May  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  0.81  cm  vertically  polarized  wavelength. 
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Figure  A-3.  Areal  map  of  the  10  May  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  1.36  cm  horizontally  polarized  wavelength. 
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Figure  A-4.  Areal  map  of  the  10  May  1979  5MMR  brightness  temperatures  (K)  for 
the  study  area  for  the  1.36  cm  vertically  polarized  wavelength. 
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Areal  map  of  the  10  May  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  1.66  cm  horizontally  polarized  wavelength. 
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Figure  A-6.  Areal  map  of  the  10  Hay  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  1.66  cm  vertically  polarized  wavelength. 
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Figure  A-8.  Areal  map  of  the  10  May  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  2.80  cm  vertically  polarized  wavelength. 
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Figure  A-9.  Areal  map  of  the  10  May  1979  SMMR  brightness  temperatures  (K)  for 
the  study  area  for  the  4.54  an  horizontally  polarized  wavelength. 
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Figure  A- 1 1 .  Areal  map  of  the  10  May  1979  maximun  temperatures  (C)  for  the  study  area 
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CORRELATION  COEFFICIENT  (xlOOO)  BETWEEN  ALL  TRANSFORMS  AND 
API  FOR  ROW  18  COL  27  BETWEEN  25  OCT  1978  AND  10  NOV  1979 
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CORRELATION  COEFFICIENT  (xIOOO)  BETWEEN  ALL  TRANSFORMS  AND 
API  FOR  ROW  21  COL  14  BETWEEN  25  OCT  1978  AND  10  NOV  1979 
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4.54  Th/T  -875  -557  -794  -826  -630  -706  -717  -273  -874  -657 

Tv-Th  743  460  538  391  584  836  561  576  872  608 

[(Tv-Th)/0.5(Tv+Th)]  766  501  594  532  602  830  722  630  907  622 

M.L.  Reg  on  All  3  899  671  863  979  702  862  922  697  952  661 


4.54  Th/T  -634  -558  -661  -875  -746  -948  -908  -882  -783  -708 

Tv-Th  791  573  628  951  344  867  554  390  249  524 

[(Tv-Th)/0.5(Tv+Th)]  777  579  687  948  432  881  545  500  387  600 

M.L.  Reg  on  All  3  843  594  902  954  924  955  921  969  934  755 
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TABLE  C-l 

BRIGHTNESS  TEMPERATURE  CORRELATION  COEFFICIENT  MATRIX 
FOR  ROW  15  COL  29  8ETWEEN  25  OCT  1978  AND  10  NOV  1979 


HORIZONTAL  POLARIZATION 
WAVELENGTH  (cm) 


0.81 

1.36 

1.66 

2.80 

4.54 

0.81 

1.00 

0.99 

0.97 

0.92 

0.80 

1.36 

1.00 

0.99 

0.95 

0.84 

1.66 

1.00 

0.98 

0.87 

2.80 

1.00 

0.92 

4.54 

1.00 

VERTICAL  POLARIZATION 
WAVELENGTH  (cm) 


0.81 

1.36 

1.66 

0.81 

1.00 

1.36 

0.98 

1.00 

1.66 

0.97 

0.99 

1.00 

2.80 

0.94 

0.97 

0.99 

4.54 
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TABLE  C-2 

NORMALIZED  BRIGHTNESS  TEMPERATURE  CORRELATION  COEFFICIENT 
MATRIX  FOR  ROW  15  COL  29  BETWEEN  25  OCT  1978  AND  10  NOV  1979 


HORIZONTAL  POLARIZATION 
WAVELENGTH  (cm) 


0.81 

1.36 

1.66 

2.80 

4.54 

0.81 

1.00 

0.95 

0.93 

0.81 

0.59 

1.36 

1.00 

0.97 

0.89 

0.69 

1.66 

1.00 

0.95 

0.76 

2.80 

1.00 

0.86 

4.54 

1.00 

VERTICAL  POLARIZATION 
WAVELENGTH  (cm) 


0.81 

1.36 

1.66 

2.80 

4.54 

0.81 

1.00 

1.36 

0.87 

1.00 

1.66 

0.87 

0.95 

1.00 

2.80 

0.81 

0.92 

0.96 

1.00 

4.54 

0.75 

0.90 

0.92 

0.96 

1.00 
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TABLE  C-3 

CORRELATION  COEFFICIENT  (xlOOO)  BETWEEN  THE  HORIZONTALLY 
POLARIZED  NORMALIZED  BRIGHTNESS  TEMPERATURE  AND  API 
BETWEEN  25  OCT  1978  AND  10  NOV  1979 


GRID  CELL 

DATA 

CHANNEL 

0.81H 

1.35H 

1.66H 

2.80H 

4.54H 

R03  C27 

Night/Day 

-768 

-771 

-816 

-800 

-794 

Night 

-715 

-802 

-811 

-843 

-875 

Day 

-848 

-789 

-840 

-821 

-775 

R09  C21 

Night/Day 

-809 

-830 

-861 

-861 

-840 

Ni  ght 

-819 

-828 

-863 

-881 

-895 

Day 

-847 

-876 

-885 

-901 

-834 

R09  C29 

Night/Day 

-748 

-782 

-800 

-780 

-747 

Ni  ght 

-701 

-796 

-787 

-785 

-758 

Day 

-831 

-799 

-830 

-826 

-790 

R15  C29 

Night/Day 

-749 

-764 

-758 

-763 

-734 

Night 

-701 

-756 

-717 

-749 

-766 

Day 

-806 

-787 

-802 

-816 

-742 

R18  C27 

Night/Day 

-678 

-719 

-736 

-736 

-752 

Ni  ght 

-600 

-707 

-685 

-685 

-734 

Day 

-771 

-760 

-787 

-807 

-787 

R21  C14 

Night/Day 

-523 

-479 

-613 

-616 

-661 

Ni  ght 

-586 

-552 

-676 

-697 

-761 

Day 

-534 

-495 

-615 

-639 

-670 

R29  C22 

Night/Day 

-647 

-660 

-693 

-708 

-720 

Night 

-658 

-678 

-687 

-763 

-747 

Day 

-679 

-672 

-716 

-712 

-758 

! 


I 
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TABLE  C-4 

CORRELATION  COEFFICIENT  (xlOOO)  BETWEEN  THE  VERTICALLY 
POLARIZED  NORMALIZED  BRIGHTNESS  TEMPERATURE  AND  API 
BETWEEN  25  OCT  1978  AND  10  NOV  1979 


GRID  CELL  DATA 


0.81V 


R03 

C27 

Night/Day 

-418 

Night 

-212 

Day 

-686 

R09 

C21 

Night/Day 

-534 

Night 

-487 

Day 

-709 

R09 

C29 

Night/ Day 

-384 

Night 

-295 

Day 

-641 

R15 

C29 

Night/Day 

-499 

Night 

-526 

Day 

-655 

R18 

C27 

Night/Day 

-518 

Night 

-462 

Day 

-692 

R21 

Cl  4 

Night/Day 

-410 

Night 

-554 

Day 

-454 

R29 

C22 

Ni ght/Day 

-359 

Night 

-458 

Day 

-414 

CHANNEL 


1.35V 

1.66  V 

2.80V 

4.54V 

-585 

-658 

-696 

-723 

-500 

-576 

-689 

-781 

-759 

-804 

-829 

-794 

-594 

-713 

-698 

-733 

-539 

-679 

-709 

-804 

-751 

-834 

-823 

-739 

-540 

-617 

-652 

-681 

-501 

-546 

-651 

-698 

-733 

-788 

-812 

-779 

-579 

-674 

-659 

-656 

-548 

-700 

-748 

-767 

-758 

-762 

-753 

-735 

-584 

-613 

-621 

-652 

-534 

-544 

-566 

-660 

-734 

-739 

-756 

-724 

-505 

-505 

-516 

-575 

-745 

-676 

-709 

-754 

-494 

-525 

-560 

-587 

-446 

-519 

-562 

-479 

-582 

-639 

-708 

-581 

-456 

-497 

-568 

-646 
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TABLE  C-5 

CORRELATION  COEFFICIENT  (xlOOO)  BETWEEN  THE  1.66  CM 
TRANSFORMS  AND  API  FOR  A  JAN-FEB  SERIES  BETWEEN 
25  OCT  1978  AND  10  NOV  1979 


GRID 

mr  TRANSFORM  MONTHS 


MAR- 

MAY- 

JUL- 

SEP- 

NOV- 

ALL 

APR 

JUN 

AUG 

OCT 

DEC 

R03 

C27 

Th/T 

-844 

-892 

-905 

-806 

-848 

-733 

Tv-Th 

783 

856 

803 

605 

735 

735 

[(Tv-Th)/0.5(Tv+Th)] 

806 

881 

795 

616 

786 

759 

M.L.  Reg  on  All  3 

873 

919 

913 

806 

953 

885 

R09 

C21 

Th/T 

-872 

-912 

-722 

-829 

-919 

-499 

Tv-Th 

813 

823 

615 

835 

886 

538 

[(Tv-Th)/0.5(Tv+Th)] 

828 

824 

649 

837 

910 

570 

M.L.  Reg  on  A1 1  3 

888 

921 

740 

853 

979 

718 

R09 

C29 

Th/T 

-814 

-956 

-347 

-812 

-819 

-849 

Tv-Th 

750 

930 

760 

718 

697 

802 

[ (Tv-Th)/0.5(Tv+Th)] 

763 

931 

778 

734 

728 

819 

M.L.  Reg  on  All  3 

820 

957 

810 

858 

859 

888 

R15 

C29 

Th/T 

-761 

-915 

-504 

-526 

-926 

-659 

Tv-Th 

679 

851 

174 

381 

872 

584 

[(Tv-Th)/0.5(Tv+Th)] 

702 

866 

220 

380 

898 

610 

M.L.  Reg  on  A1 1  3 

782 

925 

607 

529 

956 

773 

R18 

C27 

Th/T 

-734 

-921 

-461 

-753 

-946 

-692 

Tv-Th 

669 

854 

315 

730 

788 

680 

[ (Tv-Th)/0.5(Tv+Th)  ] 

687 

859 

387 

732 

827 

694 

M.L.  Reg  on  All  3 

741 

927 

781 

755 

989 

740 

R21 

C14 

Th/T 

-633 

-708 

-812 

-649 

-597 

-825 

Tv-Th 

584 

635 

652 

681 

500 

881 

[(Tv-Th)/0.5(Tv+Th)] 

587 

657 

676 

673 

505 

892 

M.L.  Reg  on  A1 1  3 

638 

748 

821 

698 

630 

907 

o* 

CM 

ex 

C22 

Th/T 

-704 

-633 

-704 

-789 

-700 

-667 

Tv-Th 

620 

612 

571 

466 

647 

713 

C(Tv-Th)/0.5(Tv+Th)] 

635 

606 

631 

497 

663 

739 

M.L.  Reg  on  All  3 

713 

681 

861 

842 

757 

795 
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TABLE  C-6 

CORRELATION  COEFFICIENT  (xlOOO)  BETWEEN  THE  1.66  CM 
TRANSFORMS  AND  API  FOR  A  FEB-MAR  SERIES  BETWEEN 
25  OCT  1978  AND  10  NOV  1979 


GRID 

CELL  TRANSFORM  MONTHS 


FEB- 

APR- 

JUN- 

AUG- 

OCT- 

ALL 

MAR 

MAY 

JUL 

SEP 

NOV 

R03 

C27 

Th/T 

-844 

-813 

-734 

-728 

-800 

Tv-Th 

783 

704 

631 

376 

735 

[ (Tv-Th )/0. 5 (Tv+Th)] 

806 

731 

651 

381 

770 

M.L.  Reg  on  All  3 

873 

849 

765 

730 

891 

R09 

C21 

Th/T 

-872 

-896 

-790 

-640 

-868 

Tv-Th 

813 

836 

811 

297 

802 

[(Tv-Th)/0.5(Tv+Th)] 

828 

841 

819 

300 

832 

M.L.  Reg  on  A1 1  3 

888 

897 

841 

643 

951 

R09 

C29 

Th/T 

-814 

-220 

-823 

-689 

-808 

Tv-Th 

750 

374 

810 

547 

674 

[(Tv-Th)/0.5(Tv+Th)] 

763 

386 

822 

548 

703 

M.L.  Reg  on  All  3 

820 

431 

866 

781 

854 

R15 

C29 

Th/T 

-761 

-744 

-442 

-422 

-799 

-836 

Tv-Th 

679 

728 

263 

208 

755 

739 

[(Tv-Th)/0.5(Tv+Th)] 

702 

753 

290 

202 

757 

762 

M.L.  Reg  on  A1 1  3 

782 

787 

502 

463 

811 

878 

R18 

C27 

Th/T 

-734 

-666 

-558 

-804 

-939 

-855 

Tv-Th 

669 

611 

475 

767 

894 

715 

c (Tv-Th ) /O . 5 ( Tv+Th ) 3 

687 

622 

503 

774 

902 

749 

M.L.  Reg  on  All  3 

741 

671 

634 

806 

942 

872 

R21 

C14 

Th/T 

-633 

-649 

-728 

-683 

-651 

-842 

Tv-Th 

584 

558 

456 

538 

658 

856 

[ (Tv-Th )/0. 5 (Tv+Th)] 

587 

578 

480 

565 

637 

866 

M.L.  Reg  on  A1 1  3 

638 

681 

791 

781 

708 

890 

R29 

C22 

Th/T 

-704 

-694 

-686 

-796 

-887 

-732 

Tv-Th 

620 

853 

479 

487 

708 

755 

[  (Tv-Th  )/0. 5  (Tv+Th)  ] 

635 

845 

548 

523 

722 

762 

M.L.  Reg  on  All  3 

713 

854 

792 

796 

891 

780 
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